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ION  BEAM  PROCESSING 


I.  INTRODUCTION 

The  primary  purpose  of  this  memo  is  to  introduce  ion  implantation  processing  to  materials  scien¬ 
tists  and  engineers  who  have  no  prior  experience  with  surface  modification  with  ion  beams.  The  intro¬ 
duction  and  background  sections  inform  the  reader  as  to  the  kinds  of  applications  that  might  be  con¬ 
sidered  for  ion  implantation  processing  and  give  a  brief  description  of  the  apparatus  and  the  physical 
principles  underlying  the  techniques.  In  subsequent  sections,  the  entire  process  of  implantation  pro¬ 
cessing  is  described  in  enough  detail  that  the  reader  will  have  a  sufficient  grasp  of  the  technology  to  be 
able  to  plan  treatment  parameters  for  parts  and/or  to  double-check  treatment  parameters  recommended 
by  vendors.  Thus,  this  memo  explains  what  ion  implantation  technology  is,  how  to  actually  implant 
ions  into  surfaces  and  gives  examples  of  parts  where  it  has  been  successfully  used. 

A  secondary  purpose  of  the  chapter  is  to  provide  workers  already  in  the  field  with  a  handbook  of 
tables,  graphs  and  equations  that  are  needed  for  doing  ion  implantation  processing.  A  few  of  the 
expressions  or  tables  are  original  but  most  of  the  information  contained  in  the  chapter  has  appeared 
elsewhere  in  the  literature  and  is  collected  here  for  convenient  reference. 

Ion  implantation  processing  of  metals  is  now  a  viable  alternative  treatment  for  specific  wear  and 
corrosion  applications,  with  vendors  distributed  in  the  U.S.,  U.K.,  Europe  and  Japan  that  provide  ser¬ 
vice  A’here  parts  can  be  treated.  The  gross  dollar  amount  of  business  for  implantation  processing  of 
metals  is  negligible  when  compared  to  other  treatments  such  as  plasma  deposition,  thermal  spraying, 
sputter  deposition,  etc.,  but  the  field  has  generated  great  interest  and  is  developing  along  a  typical 
growth  curve  for  a  high  technology  industry.  The  primary  function  of  the  vendors  is  to  treat  surfaces 
for  companies  that  are  doing  research  and  development  of  a  proprietary  nature  that  involves  the 
implantation  of  a  limited  number  of  parts  for  industrial  trials.  Some  applications  have  reached  the  pro¬ 
duction  stage  where  a  large  number  of  parts  have  been  treated  and  it  is  anticipated  that  the  number  of 
applications  reaching  production  will  increase. 

Why  is  ion  implantation  technology  useful?  One  way  to  answer  this  question  is  by  viewing  the 
entire  field  with  the  aid  of  Fig.  (1.1).  It  presents  a  periodic  table  of  the  elements  which  is  coded  with 
letters  in  the  corners  to  indicate  if  the  implanted  element  is  used  to  improve  oxidation  (O),  aqueous 
corrosion  (C),  wear  (W)  or  fatigue  (F).  In  Fig.  (1.1),  all  metals  that  have  been  implanted  are  included 
as  the  substrate  base  (e.g.,  iron  and  ferrous  alloys,  Ti  and  Ti  alloys,  etc.)  and  the  ions  cited  in  the  fig¬ 
ure  are  only  those  which  had  the  greatest  effect.  Several  features  of  this  periodic  chart  are  worth  not¬ 
ing:  i)  some  elements  improve  more  than  one  property,  ii)  Ion  implantation  is  a  versatile  process 
because  it  can  forcibly  inject  any  element  in  the  periodic  table  into  a  surface,  iii)  Elements  from  nearly 
all  the  groups  of  the  periodic  table  can  have  beneficial  effects  on  a  given  property,  iv)  Ions  which 
improve  properties  are  highlighted  here,  but  ions  which  have  deleterious  effects  may  also  be  implanted 
which  facilitates  the  study  of  mechanisms  of  wear  and  corrosion,  v)  Elements  to  implant  may  be 
chosen  for  a  particular  property  with  respect  to  the  substrate  such  as  atomic  size,  electronegativity, 
whether  the  oxides  are  n  or  p  type  semiconductors  or  insulators,  catalytic  tendency,  to  name  a  few,  and 
in  most  cases  the  element  will  form  a  solid  solution  alloy  independent  of  the  rules  of  thermodynamic 
equilibrium.  Thus,  the  interest  in  ion  implantation  lies  in  the  unparalleled  freedom  of  choice  in  the 
selection  of  ions  to  treat  and/or  study  various  surface  sensitive  properties. 

This  memo  first  informs  the  reader  of  the  types  of  applications  under  investigation  or  in  use  and 
then  introduces  the  general  features  of  ion  implantation.  This  includes  a  brief  description  of  the  pro¬ 
cess  of  injecting  ions  into  the  surface,  a  description  of  implantation  equipment,  comparisons  with  other 
surface  treatments,  and  a  complete  accounting  of  the  advantages  and  disadvantages  of  ion  implantation 
processing  Section  III  discusses  the  physical  processes  of  radiation  damage  and  phase  formation  that 
occur  during  high  dose  implantation  and  provides  a  simple  picture  of  the  development  of  microstruc¬ 
tures  that  enhance  mechanical  properties  Section  IV  gives  a  detailed  account  of  the  procedures  used  to 
ion  implant  metal  parts  including:  calculation  of  the  concentration  versus  depth  of  the  implanted  ions; 
manipulation  and  cooling  of  the  workpiece,  and  estimation  of  the  time  and  cost  of  the  treatment.  The 
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last  section  provides  case  studies  where  implantation  treatments  are  planned  and  the  treatment  parame¬ 
ters  selected  in  a  step-by-step  procedure  that  illustrates  the  use  of  the  formulas  and  tables  that  are  given 
in  the  chapter. 

II.  ION  IMPLANTATION  BACKGROUND 

A.  Applications 

Ion  implantation  is  used  in  applications  and  research  to  modify  many  surface  sensitive  properties. 
Table  (2.1)  is  a  partial  list  of  these  properties  and  shows  that  ion  implantation  can  be  applied  to  many 
materials  problems.  Although  metals  are  emphasized  in  this  chapter,  ion  implantation  is  used  to 
modify  other  materials  such  as  semiconductors,  ceramics,  insulators  and  polymers  because  it  is  univer¬ 
sally  applicable  to  any  solid  material. 

Table  2.1  —  Surface  Sensitive  Properties  Modified  by  Ion  Implantation 

Wear  Corrosion  Resistivity 

Friction  Oxidation  Dielectric  Constant 

Fatigue  Catalysis  Superconductivity 

Hardness 

Adhesion 

Toughness 


A  partial  list  of  applications  of  ion  implantation  to  metals  is  given  in  Table  (2.2).  From  this  list 
one  can  acquire  a  good  overall  view  of  the  types  of  parts  that  may  be  treated  using  ion  implantation  and 
the  kinds  of  problems  to  which  it  is  being  applied.  Some  general  points  of  interest  are  worth  noting,  i) 
the  workpieces  that  are  treated  with  ion  implantation  are  not  large,  seldom  exceeding  10  inches  in 
length  or  diameter,  ii)  the  workpieces  usually  have  a  polished  surface  finish,  iii)  most  applications  are 
for  improvement  of  wear  or  corrosion  resistance,  iv)  many  items  are  high  technology,  expensive  parts 
or  tooling.  These  points  show  that  the  focus  of  applications  is  on  small,  high-cost,  critical  parts  or  tools 
used  in  the  manufacturing  of  parts. 

Table  2.2  —  Selected  Examples  of  Applications  of  Ion  Implantation 


Paper  slitting  knives 
Rubber  slitting  knives 
Taps  for  plastics 
Extrusion  spinnerettes 
Extrusion  nozzles 
Extrusion  dies 

Pellet  punches  for  nuclear  fuel 
Forming  dies  for  cans 
Scoring  dies  for  cans 

Rolling  elements  bearings  for  air  and  space 
Instrument  bearings 
Gas  bearings 

Wire  drawing  dies  for  Cu,  Mo 
Dental  drills 


Printed  circuit  board  drills 
Wire  guides 

Orthopedic  implants  (hip  and  knee  joints,  etc.) 
Finishing  rolls  for  copper  rod 
Punches  for  acetate  sheet 
Tool  inserts 

Fuel  injectors  and  metering  pumps 
Cam  followers 
Cutters  of  plastics 

Precision  punches  for  electronic  parts 
Swaging  dies  and  press  tools  for  wheels 
Injection  mold  nozzles  and  screws 
Gears  for  ultra  high  vacuum  manipulators 


Every  surface  modification  technique  has  advantages  and  disadvantages  that  ultimately  determine 
its  applicability.  Table  (2.3)  is  a  list  of  criteria  for  when  ion  implantation  should  be  considered  as  the 
treatment  of  choice,  and  it  is  based  upon  the  advantages  and  disadvantage  that  are  discussed  in  detail 
later  in  this  section.  If  any  one  of  these  criteria  or  features  are  important  then  ion  implantation  bears  a 
close  look  as  a  possible  solution,  and  further  reading  of  this  chapter  is  recommended. 


Table  2.3  —  Criteria/Features  Influencing  the  Applicability 
of  Ion  Implantation 

Ion  implantation  should  be  considered  when  - 

1 .  Low  process  temperature  is  required  (although  not  a  necessary  condition) 

2.  Maintenance  of  high  tolerances  and/or  surface  finish  is  desired 

3.  Superior  adherence  of  treated  surface  is  important 

4.  Reproducibility  and  control  of  process  is  necessary 

5.  Parts  to  be  treated  are  small,  high-cost,  and  critical  to  operation  of  a  system 

6.  Surface  layer  of  0.1  to  1.0  microns  remains  intact  during  service  life  (note  that  lifetime  is 
frequently  improved  by  factors  of  2  to  100  by  treatment) 

7.  There  is  a  line-of-sight  access  to  the  surfaces  requiring  treatment  (a  necessary  condition) 

B.  Description  of  Ion  Implantation  Process 

Ion  implantation  is  a  process  by  which  virtually  any  element  can  be  injected  into  the  near-surface 
region  of  any  solid  by  causing  a  beam  of  high-velocity  ions,  usually  10  to  400  keV  in  energy,  to  strike  a 
target  mounted  in  a  vacuum  chamber.  The  incident  ions  come  to  a  stop  at  depths  of  tens  to  thousands 
of  angstroms  (i.e.,  0.01  to  1  /un)  in  the  host  material  as  a  result  of  losing  energy  during  collisions  with 
substrate  atoms.  The  resulting  depth  concentration  profile  of  implanted  dopant  atoms  can  be  calculated 
for  any  projectile-target  combination  from  well  established  theoretical  considerations.  At  low  ion  doses 
(i.e.,  the  number  of  ions  per  unit  area),  the  depth  concentration  profiles  are  well  characterized  by  a 
Gaussian  distribution  centered  about  an  average  projected  range.  During  the  slowing  down  process, 
collisions  between  the  incident  projectile  ions  and  the  target  atoms  result  in  the  displacement  of  target 
atoms.  This  process  is  called  radiation  damage  and  will  be  discussed  later  in  more  detail.  At  high  flu- 
ences,  where  the  concentration  of  implanted  atoms  approaches  tens  of  atomic  percent,  other  effects 
such  as  sputtering  and  ion  beam  induced  migration  of  atoms  can  significantly  alter  or  limit  the  ultimate 
concentrations  attainable.  There  is  a  probability  of  atoms  being  ejected  (sputtered)  from  the  surface  of 
the  target  as  a  result  of  collisions  in  the  surface  region,  especially  for  heavier  mass  incident  projectiles. 
An  equilibrium  condition  is  eventually  reached  where  as  many  implanted  atoms  are  removed  by 
sputtering  as  are  replenished  by  implantation.  The  depth  distribution  of  implanted  atoms  under  this 
condition  typically  has  a  maximum  at  the  surface  and  falls  off  over  a  distance  comparable  to  the  initial 
ion  range.  Figure  (2.1)  is  a  schematic  demonstration  of  the  development  of  the  atomic  concentration 
profile  during  a  high  dose  implantation.  Depth  concentrations,  radiation  damage,  sputtering  and  phase 
formation  are  discussed  in  greater  detail  in  following  sections. 

C.  Comparison  to  Coating  Techniques 

Because  ion  implantation  is  a  process  that  modifies  surface  properties  it  is  often  compared  with 
conventional  techniques,  such  as  physical  vapor  deposition,  chemical  vapor  deposition,  and  ion  plating, 
which  are  all  used  for  enhancing  surface  properties  such  as  corrosion  or  wear  resistance.  It  shares  a 
number  of  advantages  with  these  techniques  but  also  has  basic  limitations.  Therefore,  it  is  important  to 
understand  the  differences  between  ion  implantation  and  coating  techniques,  especially  ion  plating.  The 
result  of  ion  implantation  into  materials  is  the  formation  of  a  surface  alloy  of  graded  composition  that 
possesses  no  well-defined  interface  with  respect  to  the  substrate,  in  contrast  to  a  deposited  layer  where 
the  interface  is  very  abrupt.  The  thickness  of  the  implanted  region  is  typically  less  than  2000 


Angstroms  (.3)  for  implantation  energies  of  100  keV  (1000  A  for  heavy  ions).  In  ion  plating,  how¬ 
ever,  the  coating  is  typically  much  thicker,  and  its  composition  is  independent  of  the  nature  of  the  sub¬ 
strate.  Although  ion  plating  is  often  carried  out  with  the  substrates  electrically  biased  at  several  keV, 
the  mean  energy  of  the  charged  particles  reaching  the  surface  is  of  the  order  of  only  100  eV.  This 
occurs  because  only  a  small  fraction  (—10%)  of  the  particles  are  ionized  when  they  reach  the  substrate. 
Hence,  there  is  virtually  no  penetration  of  the  incident  beam  into  the  substrate.  Ion  plating  is  carried 
out  under  a  relatively  high  pressure  that  causes  atomic  collisions  and  scattering.  This  factor  is  responsi¬ 
ble  for  the  high  "throwing  power”  of  the  ion  plating  technique  versus  the  "line-of-sight"  limitation  of 
implantation. 

D.  Historical  Perspective 

Since  the  early  1970’s,  ion  implantation  has  been  extensively  used  by  the  semiconductor  device 
industry  as  a  method  of  introducing  controlled  amounts  of  dopants  into  the  surface  region  of  semicon¬ 
ductor  substrates.  It’s  principal  advantages  for  this  application  include  improved  controllability  and 
reproducibility  for  device  fabrication,  as  compared  to  thermal  diffusion.  The  absolute  concentration  of 
implanted  atoms  and  the  uniformity  across  the  sample  surface  may  be  controlled  to  better  than  5%  and 
1%  respectively.  The  volume  concentrations  of  impurities  required  for  semiconductor  applications  typi¬ 
cally  range  from  parts  per  million  to  as  high  as  0.1  at  .%,  whereas  the  concentrations  required  for  many 
of  the  non-semiconductor  applications,  such  as  for  bearings,  are  typically  ten  to  thirty  atomic  percent. 
Development  of  techniques  and  equipment  for  ion  implantation  is  driven  by  the  semiconductor 
industry;  the  equipment  for  implantation  of  metals  is  a  spin-off  from  that  effort.  Therefore,  the  tech¬ 
nology  ion  implantation  of  metals  advances  along  with  technological  progress  in  the  semiconductor 
field. 

E.  Ion  Implantation  -  Advantages  and  Limitations  of  the  Technique 

Some  of  the  advantages  and  limitations  of  ion  implantation  in  comparison  to  other  surface  treat¬ 
ments  (such  as  coatings)  are  listed  in  Table  (2.4). 


Table  2,4  —  Advantages  and  Limitations  of  Ion  Implantation 
as  a  Surface  Modification  Technique 


Advantages 

Limitations 

(1) 

Low  temperature  process 

(2) 

No  significant  dimensional  changes 

(3) 

No  adhesion  problems  since  there 
is  no  sharp  interface 

(4) 

No  sacrifice  of  bulk  properties 

(5) 

Solid  solubility  limit  can  be 

(1) 

Line-of-sight  process 

exceeded 

(2) 

Shallow  penetration 

(6) 

Alloy  preparation  independent  of 

(3) 

Relatively  expensive 

diffusion  constants 

equipment  and  processing 

(7) 

Allows  fast  screening  of  the  effects 
of  changes  in  alloy  composition 

costs 

(8) 

Depth  concentration  distribution 
controllable 

(9) 

Highy  controllable  and  reproducible 

The  chief  advantages  of  ion  implantation  are  low  temperature  processing  (no  distortion  of  part), 
no  change  in  dimensions  (no  refinishing  required),  excellent  adhesion  of  the  implanted  surface,  and 
that  the  surface  properties  can  be  optimized  independently  from  the  bulk  properties.  The  fact  that  ion 
implantation  is  a  nonequilibrium,  fast  quench  technique  permits  the  formation  of  surface  alloys 
independently  of  solubility  limits  and  diffusi vities  that  govern  conventional  alloy  formation. 


Ion  implantation  often  allows  the  convenient  production  and  subsequent  study  of  surface  alloys 
with  well  defined  compositions.  In  this  manner  the  technique  can  be  used  as  a  powerful  research  tool 
to  examine  the  physical  state  of  alloys  as  a  function  of  varying  alloy  composition.  Using  implantation 
alloying  one  can  avoid  changing  other  parameters  (e  g.,  grain  size)  which  may  effect  the  property  of 
interest  such  as  oxidation  behavior. 

A  basic  limitation  of  ion  implantation  is  that  as  a  tine-of-sight  process,  it  is  not  feasible  to  apply  it 
to  samples  having  complicated  re-entrant  surfaces.  The  shallow  depth  of  penetration  would  seemingly 
limit  the  technique  for  engineering  applications;  however,  there  are  many  situations  involving  both 
physical  and  chemical  properties  in  which  the  effect  of  the  implanted  ions  is  to  change  the  mode  of 
wear  or  corrosion  behavior  resulting  in  significant  extension  of  lifetime  of  the  part.  The  disadvantages 
determine  the  applicability  of  a  process.  The  limitation  of  shallow  depth  and  relatively  high  cost  imply 
that  the  most  appropriate  applications  are  for  high  cost,  high  technology  critical  parts  which  limit  the 
performance  or  lifetime  of  end  use  hardware  or  expensive  tools  used  in  the  production  of  hardware. 

F.  Ion  Implantation  Equipment 

A  typical  production-type  ion  implanter  is  shown  in  Fig.  (2.2).  It  consists  of  an  ion  source  capa¬ 
ble  of  producing  ion  beams  of  practically  any  stable  element.  Ions  are  extracted  from  the  ion  source  by 
an  electrode  held  at  high  potential.  Since  ion  sources  normally  produce  ions  of  several  species  in  addi¬ 
tion  to  the  one  desired  for  implantation,  it  is  ordinarily  essential  to  mass-analyze  the  beam  in  order  to 
allow  only  the  species  of  interest  to  continue  toward  the  target.  The  separation  of  the  isotopes  of  Cr  is 
depicted  here.  After  mass  analysis,  the  beam  is  accelerated  through  an  evacuated  acceleration  column 
during  passage  from  the  high  voltage  terminal  to  ground  potential.  After  acceleration  the  beam  is  elec¬ 
trostatically  focussed  and  raster  scanned  over  an  aperture  held  in  front  of  the  target  to  ensure  unifor¬ 
mity.  The  ion  beam  current  striking  the  target  is  then  current  integrated  to  obtain  the  impurity  concen¬ 
tration  (in  dimensions  of  atoms  per  unit  area)  which  in  turn  can  be  converted  into  a  volume  concentra¬ 
tion  by  knowledge  of  the  range-energy  relationships  for  the  projectile-ion  host  combination.  A 
schematic  diagram  of  a  typical  ion  source  for  metals  implantation  is  shown  in  Fig.  (2.3).  A  low  pres¬ 
sure  plasma  is  sustained  by  an  arc  discharge  voltage  of  100V  and  by  a  hot  filament  to  supply  electrons. 
The  positive  ions  that  are  extracted  from  the  plasma  through  a  25kV  potential  can  be  composed  of  the 
gas  used  to  sustain  the  arc  discharge,  or  of  metallic  ions  derived  from  volatile  metal  chlorides  produced 
in  a  reaction  chamber  in  the  rear  of  the  ion  source.  Ion  beams  of  most  of  the  elements  in  the  periodic 
table  can  be  produced  in  this  simple  ion  source  configuration. 

Non-mass  analyzed  implanters  are  available  for  applications  requiring  ions  formed  from  pure  gases 
and  are  cheaper  to  manufacture.  The  newest  generation  commercial  implanters  are  capable  of  beam 
currents  of  10  mA,  and  development  of  100  mA  machines  is  in  progress.  The  flux  of  ions  for  these 
beam  currents  is  6.2  x  10lh  and  6.2  x  10i7  ions/s,  respectively.  A  typical  dose  for  metals  implantation 
is  2  x  1017  ions/cnr.  The  time  required  to  treat  one  square  inch  of  flat  surface  with  these  beam 
currents  would  be  21  seconds  and  2.1  seconds,  respectively. 

This  introduction  to  implantation  has  been  very  brief  because  the  process  itself  is  straight-forward 
and  therefore  easy  to  understand.  For  further  reading  on  the  equipment  or  other  aspects  of  ion  implan¬ 
tation  covered  in  this  report,  the  references  cited  in  Table  (2.5)  are  recommended. 


If 


f 

k 


Table  2.5  —  References  for  Information  Related  to  Ion  Implantation 

1 .  Applications 

P.  Sioshansi,  Thin  Solid  Films  118  (1984)61 
R.E.  Hoisington,  D.E.  Technology,  (Jan.  1984)10. 

2.  Catalysis 

G.K.  Wolf,  Nucl.  Instrum,  and  Meth.  182/183  (1981)875. 

3.  Corrosion 

C.R.  Clayton,  Nucl.  Instrum,  and  Meth.  182/183  (1981)865. 

4.  Deformation 

D  A.  Hills  and  D.W.  Ashelby,  Wear  75  (1982)221. 

5.  Equipment  and  Techniques 

G.  Dearnaley  and  P.D.  Goode,  Nucl.  Instrum,  and  Meth.  189  (1981)117. 

J. W.  Butler,  Encyclopedia  of  Chemical  Technology  (John  Wiley,  1981) 
Vol.  13,  p.  706. 

6.  Fatigue  (Ion  Implantation) 

R. G.  Vardiman,  in  "The  Application  of  Ion  Plating  and  Ion 
Implantation  to  Materials"  Ed.,  R.W.  Hockman,  ASM,  (1986). 

7.  General  Ion  Implantation  Review 

S. T.  Picraux,  Ann.  Rev.  Mater.  Sci.  14  (1984)335. 

8  Oxidation 

K. S.  Grabowdki  and  L.E.  Rehn,  in  "Corrosion  of  Metals  Processed  by 
Directed  Energy  Beams"  (eds.,  C.R.  Clayton  and  C.M.  Preece,  The 
Metallurgical  Society  of  AIME,  Louisville,  Kentucky,  1981)  p.  23. 

9.  Phase  Formation 

S.T.  Picraux,  in  "Site  Characterization  and  Aggregation  of  Implanted 
Atoms  in  Materials"  (eds.,  A.  Perez  and  R.  Coussement;  Plenum  NY, 
1980)  pp.  307  and  325. 
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10  Wear 

"Fundamentals  of  Friction  and  Wear  of  Materials,"  Ed.,  D.A.  Rigney, 
ASM.  Metals  Park,  OH  (1981). 

11  Wear  (Ion  Implantation) 

I.L.  Singer,  Materials  Research  Society  Symp.  Proc.  Vol.  27 
(Elsevier,  NY,  1984)  p.  585. 
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III.  RADIATION  DAMAGE  AND  PHASE  FORMATION 


A.  Introduction 

This  section  gives  a  brief  introduction  to  the  basic  processes  that  occur  on  an  atomic  scale  during 
ion  implantation  in  order  to  provide  an  understanding  of  how  the  composition  and  microstructure  of 
surfaces  are  modified  as  the  ion  dose  increases  and  impurity  atoms  are  forced  into  a  crystalline  host  to 
concentrations  of  10  to  30  atomic  percent. 


An  energetic  ion  loses  energy  in  its  passage  through  a  solid  via  electronic  and  atomic  collisions. 
As  far  as  damage  is  concerned,  electronic  collisions  do  nothing  more  than  heat  up  the  lattice  (Czjzek  et 
al.,  1968).  On  the  other  hand,  atomic  collisions  can  recoil  lattice  atoms  with  energies  well  in  excess  of 
the  minimum  energy  for  displacement,  ED  which  for  most  materials  is  ~  30  eV,  and  the  recoils  conse¬ 
quently  undergo  further  collisions  with  other  lattice  atoms  before  they  themselves  come  to  rest.  This 
multiplication  process  is  called  a  "collision  cascade"  and  for  incident  energies  relevant  to  ion  implanta¬ 
tion  (10  to  400  keV)  can  involve  thousands  of  atoms  and  result  in  hundreds  of  atoms  becoming  per¬ 
manently  displaced  from  equilibrium  lattice  sites.  Holes  which  remain  in  the  crystalline  lattice  after  dis¬ 
placement  are  called  vacancies,  and  the  displaced  atoms  which  become  squeezed  into  the  lattice  are 
called  interstitials.  These  vacancies,  interstitials  and  defect  aggregates  produced  by  atomic  collisions 
comprise  what  is  commonly  called  "radiation  damage."  Figure  (3.1)  illustrates  a  possible  sequence  of 
events  in  a  single  collision  cascade. 


The  maximum  energy  that  a  lattice  atom  can  receive  in  a  head-on  collision  is  large  and  is  given  by 

4 M]  M2 


£*max  _ 


(A/ 1  +  M2) 
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where  E0  is  the  incident  ion  energy  and  Mx  and  M2  are  the  ion  and  target  atomic  masses,  respectively. 
The  average  primary-knock-on  atom  (PKA)  energy  is  much  smaller  than  £max  because  the  random  dis¬ 
tribution  of  trajectories  in  the  solid  makes  large  impact  parameters  much  more  likely  than  small  ones. 
Therefore  average  PKA  energies  are  on  the  order  of  a  few  keV.  The  following  discussion  therefore 
focusses  on  the  radiation  damage  of  few  keV  PKA  iron  atoms  in  an  iron  target  in  order  to  illustrate  the 
nature  of  radiation  damage. 

B.  The  Structure  of  Radiation  Damage 


The  damage  caused  by  the  collision  cascade  of  one  primary  knock-on  atom  (PKA)  with  an  energy 
of  1  keV  (i.e.  one  iron  lattice  atom  suddenly  given  1  keV  of  energy)  can  be  estimated  from  the  modi¬ 
fied  theory  of  Kinchin  and  Pease  (Kinchin  and  Pease,  1955).  They  suggested  that  the  number  of 
Frenkel  defects  nD  (1  defect  =  1  interstitial  +  1  vacancy)  created  by  one  PKA  with  energy  Ep  can  be 
expressed  as 


nD  = 


0-8  EP 
2  Ed 


(3.1) 


where  E0  is  the  minimum  energy  for  a  displacement.  For  a  1000  eV  recoil  in  Fe  (ED  =  30  eV)  (Ergin- 
soy  et  al..  1964)  nD  equals  13.  In  an  important  work,  Erginsoy  et  al.  performed  a  computer  study  by 
means  of  Monte  Carlo  techniques  of  low  energy  knock-ons  in  iron  at  0nA  where  no  defect  annealing 
occurs  and  found  that  the  number  of  defects  created  could  be  well  described  by  equation  (3.1)  for  Er 
energies  of  <T  1500  eV  (Erginsoy  et  al.,  1965). 

Another  consequence  of  collision  cascades  is  that  they  can  lead  to  a  dynamic  separation  of  the 
vacancy  and  its  associated  interstitial  which  is  very  much  greater  than  that  for  random  collisions  at  the 
same  energy  Accordingly,  computer  studies  of  Erginsoy  et  al.  (Erginsoy  et  al.,  1964  and  1965)  have 


shown  that  the  displacement-spike  (damage  due  to  one  primary  knock-on)  has  a  non-uniform  spatial 
distribution  of  damage.  The  vacancy  population  is  higher  than  the  interstitial  in  the  central  part  of  the 
cascade,  but  the  reverse  is  found  in  the  annular  volume  surrounding  the  central  vacancy  rich  region. 
Beeler  et  al.  (Beeler,  1964  and  1969)  found  that  on  the  average,  the  structure  of  damage  for  a  1  keV 
primary  knock-on  is  a  vacancy  concentration  of  3  atomic  %,  with  monovacancies  (isolated  vacancies) 
comprising  about  35%  of  the  total  vacancy  population,  divacancies  20%,  and  clusters  of  three  or  more 
vacancies  45%.  About  3%  of  the  vacancies  are  contained  in  clusters  of  more  than  ten  vacancies.  More 
than  95%  of  the  interstitials  atoms  are  mono-interstitials  and  the  remainder  di-interstitials.  Following 
the  work  of  Beeler,  Doran  performed  a  similar  study  for  PKA’s  between  5  and  20  keV  and  found  that 
the  damage  consisted  mainly  of  isolated  interstitial  clusters  and  vacancy  clusters  (Doran,  1970). 

C.  Annealing  of  Damage 

In  iron  considerable  annealing  of  damage  takes  place  at  room  temperature.  Johnson  (Johnson, 
1964)  performed  computer  simulations  of  defects  in  an  iron  lattice  and  found  the  migration  energies  of 
the  defects  to  be  0.68  eV  for  a  monovacancy,  0.66  eV  for  a  divacancy,  0.33  eV  for  a  monointerstitial, 
and  0.18  eV  for  a  di-interstitial.  At  room  temperature  the  vacancies  are  considered  immobile,  while  the 
interstitials  are  mobile.  Clusters  of  three  or  more  vacancies  or  interstitials  are  considered  to  be  immo¬ 
bile.  Hence,  the  annealing  takes  place  when  interstitials  find  vacancies.  There  is  always  residual  dam¬ 
age  after  annealing  is  completed  since  some  of  the  mobile  interstitials  may  cluster  and  become  immo¬ 
bile  (higher  migration  energy)  or  travel  to  dislocations,  grain  boundaries  or  the  sample  surface  which 
act  as  interstitial  sinks.  This  leads  to  some  net  residual  vacancy  concentration.  At  elevated  tempera¬ 
tures  the  migration  of  interstitials  and  vacancies  is  faster  and  annealing  more  complete.  The  activation 
(migration)  energy  for  motion  Em  associated  with  each  defect  gives  roughly  the  mean  time  r,  for  one 
jump  using  the  formula, 

kT  =  0.026  eV  (at  300  K ) 

E 

T  ”  T°  e  ~kT  to  —  10~Us-  (3.2) 

Beeler  (1964)  studied  annealing  in  the  damage  spike  during  25  jumps  per  existing  mobile  interstitial, 
and  found  that  50%  of  all  mobile  interstitials  were  removed  within  these  25  jumps.  Since  the  possible 
vacancies  for  the  remainder  of  the  interstitials  to  combine  with  is  reduced,  the  number  of  jumps  needed 
to  completely  anneal  out  the  interstitials  is  much  larger  than  25.  Doran  (1970)  has  used  computer 
simulation  to  calculate  displacement  spike  annealing  at  room  temperature  for  a  PKA  energy  of  5  to  20 
keV  in  iron  using  the  migration  energies  of  Johnson.  Doran’s  studies  found  that  allowing  only  intersti¬ 
tials  to  be  mobile,  annealing  was  completed  after  500  jumps  per  interstitial,  and  that  about  70%  of  the 
damage  was  removed.  Many  interstitials  remained  as  immobile  clusters  of  three  or  more  interstitials. 
Upon  allowing  a  small  monovacancy  mobility,  annealing  was  complete  after  6000  interstitial  jumps  and 
1260  vacancy  jumps,  and  about  95%  of  the  damage  was  removed,  the  remainder  being  immobile 
vacancy  clusters.  This  would  be  likely  to  characterize  the  implantation  damage  at  low  doses  and  for 
implantation  at  room  temperature.  As  the  dose  increases,  the  individual  collision  cascades  overlap  and 
every  atom  in  the  implanted  region  is  eventually  displaced  and  moved  to  another  lattice  or  defect  posi¬ 
tion.  At  high  doses,  a  dynamic  balance  is  set  up  where  as  many  frenkel  pairs  are  annihilated  as  are 
created  and  saturation  of  the  radiation  damage  occurs. 

D.  Damage  at  High  Dose 

Expressions  for  the  total  amount  of  energy  deposited  in  displacement  collisions  (atomic  collisions) 
per  ion,  per  unit  length,  is  given  in  a  later  section.  In  Figure  (3.2)  Sn  is  plotted  for  150  keV  Fc 
ions  implanted  into  an  Fe  substrate.  The  depth  distribution  of  implanted  iron  ions  is  also  shown  which 
has  a  Gaussian  shape  and  is  characterized  by  the  mean  range  Rn  and  standard  deviation  of  the  range 
6  Rp.  Note  that  the  maximum  in  the  damage  distribution  occurs  at  a  shallower  depth  than  Rr.  The 
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number  of  displacements  vs  depth  can  be  obtained  with  the  modified  Kinchin  and  Pease  (1955) 


expression. 


nD  (x ) 


0.8  S„  (x) 
2  Ed 


and  the  total  number  of  displacements  in  the  implanted  volume  is 

nD  =  SJx)  dx. 


A  commonly  used  unit  of  damage  is  displacements  per  atom  (DPA).  A  unit  of  1  DP  A  means  that  on 
the  average,  every  atom  in  the  affected  volume  has  been  displaced  from  its  equilibrium  lattice  site 
once.  The  dependence  of  DPA  vs  depth  is  obtained  by  dividing  equation  (3.3).  the  number  of  dis¬ 
placed  atoms  per  unit  length,  by  the  atomic  density  N  (atoms/cm3)  or; 

0.8  S„  (x) 

DPA  (x )  =  -  </>  (3.4) 

2EdN 

where  <b  is  the  ion  dose  in  units  of  atoms/cm2.  For  ion  masses  A/,:>20  and  for  energies  of  10  to  400 
keV,  the  displacement  damage  profile  in  Figure  (3.2)  is  quite  typical  so  that  an  estimate  of  DPA  for 
these  ion/target  combinations  can  be  made  by  assuming  a  uniform  damage  profile  to  a  depth  Rp  and  by 
using  S„(x=0)  as  the  total  energy  expended  in  displacement  collisions.  Then  the  DPA  in  the 
implanted  region  is  simply 

0.8  SJx  =  o) 

DPA  -  - r ==-r: -  <t>  (3.5) 

2  EdN 

For  the  ions  and  energies  relevant  to  ion  implantation,  calculations  using  equation  (3.5)  indicate 
that  a  damage  dose  of  1  DPA  occurs  in  the  dose  range  of  1014  (for  mass  150)  to  10lh  (mass  8) 
atoms/cm2.  Typical  doses  for  implantation  of  metals  are  10|7/cm2  or  greater.  Hence  each  atom  in  the 
implanted  volume  is  displaced  tens  to  hundreds  of  times  and  the  individual  collision  cascade  volumes 
overlap.  For  self  ion  implantation  (e  g.  Fe  into  Fe),  transmission  electron  micrographs  show  a  high 
density  of  dislocation  tangles,  vacancy  clusters,  and  vacancy  and  interstitial  loops  as  contrasted  with  low 
dose  (<  <  1  DPA)  where  individual  vacancy  clusters  are  observed  that  are  produced  in  a  single  collision 
cascade.  The  excess  vacancy  concentrations  from  individual  cascades  coalesce  into  these  defect  struc¬ 
tures  as  the  dose  increases.  It  has  often  been  observed  that  the  highly  ion  damaged  layer  resembles  a 
surface  that  has  undergone  severe  cold  work,  and  by  analogy,  the  surface  should  be  in  a  state  of 
compression.  This  is  certainly  true  if  the  implanted  ion  has  a  larger  atomic  size  than  the  substrate. 

E.  Final  Site  of  Implanted  Ions 

The  computer  simulation  work  of  Erginsoy  et  al.  (1965)  predicts  that  for  low  energy  PKA's  in 
iron  the  primary  always  ends  up  in  a  substitutional  site  via  a  'eplacement  collision  with  the  associated 
interstitial  several  sites  away,  and  the  vacancy  at  a  first,  second,  or  third  nearest  neighbor  site  Dau- 
treppe  (Dautreppe,  1971)  has  reviewed  the  question  of  the  final  site  of  the  implant  and  found  that 
experiment  and  theory  suggest  the  findings  of  Erginsoy  et  al.  are  correct  as  long  as  .V/,  >  V/:  There 
are  expectations  such  as  the  cases  of  Xe  (Feldman  and  Murnick.  1971)  and  Pb  (Feldman  et  al..  1 97 1 ) 
in  iron.  In  general,  each  case  of  Mx  in  any  target  should  be  examined  individually  to  be  certain  of  the 
final  site  of  the  implant. 

There  are  a  number  of  papers  that  report  lattice-site  location  of  implanted  ions  in  metals  (Vian- 
den,  1980;  Borders  and  Poate,  1976;  Sood,  1978).  Borders  and  Poate  state  that  for  concentrations  less 
than  or  equal  to  1  at  .%,  substitutional  and  interstitial  solutions  have  been  observed  for  implantation 
into  Be,  Fe,  Ni  and  Cu,  and  at  concentrations  greater  than  10  at  both  metastable  solid  solutions 


and  amorphous  alloys  have  been  formed  in  Fe,  Ni  and  Cu.  The  substitutional  fraction  tends  to  increase 
with  increasing  atomic  number  of  implanted  ions  (Borders  and  Poate,  1976).  For  example,  it  was 
found  that  B  implanted  into  Ni  occupies  interstitial  sites  and  Sn  implanted  into  nickel  is  located  at  sub¬ 
stitutional  sites  with  a  high  probability  (>60%). 

Light  ions  such  as  H,  He,  B,  C,  N,  and  0  usually  reside  at  interstitial  sites  in  the  lattice  unless 
compound  formation  changes  the  phase  or  structure  of  the  target.  This  can  occur  at  high  doses  when 
the  implanted  ion  concentration  exceeds  a  10  at.%. 

F.  Phase  Formation 

A  very  useful  picture  to  represent  the  temporal  dependence  of  a  collision  cascade  has  recently 
been  proposed  by  Johnson  et.  al,  to  explain  the  phenomenon  of  ion  mixing  (Johnson  et  al.,  1985).  We 
use  this  model  to  introduce  the  factors  influencing  phase  formation  because  of  the  strong  similarity 
between  ion  mixing  of  bi-layer  and  multilayer  metallic  samples  and  high  dose  ion  implantation.  In  the 
latter,  ion  mixing  processes  are  occurring  when  the  impurity  concentration  becomes  greater  than  —10 
at  .%  because  subsequent  implanted  atoms  ion  mix  the  previously  implanted  atom  concentration. 

The  time  evolution  of  a  system  undergoing  ion  mixing  may  be  divided  into  two  regimes.  The 
prompt  regime  (PR)  is  the  time  from  the  initial  interaction  to  when  thermal  equilibrium  is  rees¬ 
tablished  within  the  metal  (t  =  0  to  —  10~9s).  Longer  times  are  called  the  delayed  regime  (DR)  where 
nonequilibrium  concentrations  of  defects  may  exist  which  can  continue  to  migrate  in  the  delayed 
regime  (/>  10  9s)  as  earlier  described  by  the  calculations  of  Doran  and  Beeler.  Figure  (3.3)  illustrates 
the  various  regimes.  Following  the  work  of  Johnson  et  al.  (1985),  the  prompt  regime  may  be  further 
broken  down  into  two  sub-regimes.  The  time  required  for  the  ion  and  secondary  ions  in  the  collision 
cascade  to  lose  most  of  their  energy  is  defined  as  the  ballistic  regime  and  occurs  in  the  approximate 
time  regime  from  the  initial  interaction  to  10-13  seconds.  The  time  required  for  the  kinetic  energy  to 
be  shared  equally  among  the  atoms  in  the  volume  of  the  cascade  where  a  thermal  equilibrium  is  es¬ 
tablished  is  defined  as  the  thermalizing  regime  and  occurs  in  the  approximate  time  regime  of  10-13  to 
10  9  seconds.  The  majority  of  atomic  mixing  occurs  in  the  ballistic  regime  and  in  the  initial  stages  of 
the  thermalizing  regime  where  high  local  equilibrium  temperatures  (perhaps  even  melting)  can  cause 
diffusional  mixing  over  distances  on  the  order  of  the  cascade  volume.  The  amount  of  mixing  that 
occurs  in  the  ballistic  regime,  0  to  10~13  seconds,  depends  to  First  order  upon  the  mass  and  atomic 
charge  of  the  ion,  the  energy  of  the  ion,  and  the  atomic  masses  and  atomic  charges  of  the  two  metals 
involved  Mixing  in  the  thermalizing  regime  in  the  time  span  of  10~13  to  10“"  seconds  can  be  influ¬ 
enced  by  thermodynamics  as  has  been  shown  by  Cheng  et  al.  (Cheng  et  al.,  1985).  Heats  of  mixing 
play  a  role  in  determining  the  rate  constant  for  mixing  kinetics  which  proceeds  as  ion  dose  to  the  1/2 
power 

Most  of  the  mixing  of  A  atoms  with  B  atoms  and  therefore  the  establishment  of  the  atomic  con¬ 
centrations  occurs  in  the  ballistic  regime  and  the  thermalizing  regime  from  time  0  to  10~u  seconds. 
Phase  formation  occurs  in  the  thermalizing  regime  in  the  time  span  of  10-11  to  10~9  seconds  where  the 
average  energy  of  the  atoms  in  the  cascade  goes  from  I  eV  to  thermal  equilibrium  values.  Any  phases 
which  form  must  arise  from  diffusionless  transformations  or  from  near  epitaxial  growth  on  existing  «- 
solution  or  /j-solution  surfaces.  This  rapid  quench  rate  rules  out  the  formation  of  phases  with  complex 
structures  lor  compositions  where  a-  and  /J-solulions  are  unstable  to  melting  or  amorphization,  amor¬ 
phous  phases  are  formed 

In  the  phase  formation  time  regime  of  10  11  to  10“9  seconds  nucleation  theory  of  solid-solid  or 
liquid-solid  transformations  can  he  used  to  rationalize  the  influence  of  thermodynamic  driving  forces  on 
the  phases  which  form  Figure  (3  4)  is  a  temperature-time  to  transformation  plot  illustrating  the  kinet¬ 
ics  of  a  solid-solid  phase  transformation  upon  cooling  from  the  equilibration  temperature  Teq.  At  tem¬ 
peratures  near  Trq.  the  rate  of  formation  of  the  new  phase  is  controlled  by  nucleation  and  epitaxial 
growth  As  the  temperature  drops,  diffusion  of  atoms  to  the  phase  interfaces  controls  the  rate  of  phase 


formation.  The  lower  curve  shows  that  for  a  linear  temperature  drop  with  time  at  a  slow  cooling  rate, 
the  temperature  passes  though  a  region  where  growth  of  the  phase  is  allowed  and  complete  transforma¬ 
tion  to  the  new  phase  occurs.  For  the  faster  cooling  rate,  the  region  of  formation  of  the  new  phase  is 
completely  bypassed  and  metastable  or  amorphous  phases  are  likely  to  be  produced.  For  the  phase  for¬ 
mation  region  in  Figure  (3.3)  the  temperature  drops  by  103  to  104K  in  10“ 10  to  10“ 11  seconds  so  that 
the  effective  quench  rate  is  on  the  order  of  1014  K/s.  This  is  much  faster  than  melt-quenching  tech¬ 
niques  ( 1 07 — 109  K/s).  The  primary  differences  between  quenching  in  a  cascade  and  melt  quenching  is 
that  cascade  quenching  occurs  in  the  solid  state,  is  accompanied  by  a  high  density  of  defects,  and  is  fas¬ 
ter.  Further  approach  to  equilibrium  can  occur  in  the  delayed  regime  because  the  presence  of  a  high 
density  of  defects  that  can  assist  in  atomic  transport.  Thus,  precipitation  and  segregation  that  began  in 
the  phase  formation  regime  can  continue  in  the  delayed  regime.  Formation  or  ripening  of  phases  in  the 
regime  is  controlled  by  the  mobility  of  defects  and  defect  complexes  and  therefore  differs  widely  among 
materials. 

Interstitial  and  substitutional  metastable  solid  solutions  are  frequently  observed  in  ion  implanted 
surfaces.  Molybdenum  which  has  negligible  solubility  in  Al,  can  be  formed  in  continuous  metastable 
solid  solution  in  Al  up  to  20  at.%  concentration  (Benson,  1984).  Nickel-Boron,  Fe-B  and  Fe-P  amor¬ 
phous  alloys  prepared  by  melt  quenching  are  also  readily  formed  by  implantation  of  the  metalloid  into 
the  metal  (Grant,  1981).  The  ability  to  form  unique  metastable  and  amorphous  alloys  is  one  of  the 
most  attractive  features  of  ion  implantation  processing.  Studies  cited  earlier  in  the  chapter  (Table  2.5) 
indicate  that  metastable  surface  alloys  possess  interesting  and  useful  properties  for  tribiological  and  cor¬ 
rosion  applications. 


IV.  ION  BEAM  PROCESSING 


A.  Introduction 

Once  it  is  decided  that  the  implantation  of  a  certain  ion  is  required  to  treat  a  part,  questions  arise 
such  as  what  is  the  correct  ion  dose  and  energy?  What  incident  angle  of  the  ion  beam  should  be  used 
to  implant  facets  or  curved  surfaces  on  the  workpiece?  How  is  the  workpiece  cooled  and  mounted  in 
the  vacuum  chamber?....  and  so  on.  This  section  organizes  answers  to  these  and  other  questions  into 
manageable  elements  and  is  designed  to  allow  the  reader  to  make  informed  estimates  of  the  ion  beam 
parameters  that  should  be  used  for  proper  treatment  of  a  workpiece.  These  considerations  for  choosing 
ion  beam  parameters  are  similar  to  those  used  by  ion  implantation  vendors.  By  being  familiar  with  the 
principles  and  terminology,  it  should  facilitate  discussions  with  vendors  and  allow  the  materials  engineer 
to  double-check  the  implantation  conditions  that  vendors  recommend. 

B.  Ion  Range  Parameters 

i.  General  Features  of  Range  Distributions 

The  most  frequently  asked  question  about  ion  implantation  is  "what  is  the  depth  of  the  implanted 
layer"?  This  section  discusses  the  factors  determining  the  ion  range  and  range  distribution  and  includes 
a  table  from  which  the  ion  range  can  be  conveniently  estimated  for  many  ion/target  combinations. 

The  depth  of  penetration  (ion  range)  Rp  varies  between  about  5  nm  and  1  for  the  ion/target 
combinations  of  interest,  with  50  nm  being  rather  typical.  Penetration  occurs  by  means  of  momentum 
expenditure  rather  than  by  thermal  diffusion  into  the  lattice.  Therefore,  the  ion  range  is  small  com¬ 
pared  to  the  thickness  of  diffusion  coatings  (carburization,  nitriding,  etc.).  Figure  (4.1)  is  a  plot  of  the 
calculated  projected  range  Rp  and  standard  deviation  of  the  projected  range  &Rr  (also  called  straggling) 
of  nitrogen  in  iron  versus  the  incident  nitrogen  ion  energy.  As  shown  in  the  insert,  the  range  distribu¬ 
tion  of  implanted  ions  is  well  represented  by  a  Gaussian  curve  at  low  doses.  Figure  (4.2)  shows  the 
calculated  dependence  of  the  ion  range  distribution  on  initial  ion  energy  for  iron  incident  on  an  iron 
substrate  (Butler,  1981).  The  implantation  doses  are  equal  for  all  four  cases.  As  the  energy  increases 
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the  ion  range  Rp  increases,  the  width  of  the  distribution  A/?p  increases,  and  the  peak  concentration 
decreases.  Consequently,  the  dose  for  the  200  keV  implant  would  have  to  increase  in  order  to  maintain 
the  same  peak  concentration  as  for  lower  energy  implants.  This  general  behavior  applies  to  all 
ion/ target  combinations.  Two  additional  trends  are  worth  noting.  For  a  given  substrate,  Rp  decreases 
as  the  atomic  number  of  the  ion  increases;  for  a  given  ion,  Rp  also  decreases  as  the  atomic  number  of 
the  substrate  increases.  The  latter  trend  depends  upon  the  substrate  density  and  significant  departures 
from  this  rule  are  possible  in  regions  of  the  periodic  table  where  the  density  of  the  elements  changes 
rapidly. 

A  Gaussian  shape  for  the  profile  arises  from  the  statistical  nature  of  the  stopping  process  wherein 
the  ion  loses  energy  in  collisions  with  both  electrons  and  atoms  in  the  solid.  Only  a  small  amount  of 
energy  is  exchanged  per  collision  with  electrons  because  of  the  small  electron  mass.  A  much  larger 
amount  of  energy  per  collision  is  transferred  in  collisions  with  the  more  massive  atoms  in  the  solid  but 
the  absolute  value  of  the  average  energy  per  atomic  collision  is  still  much  less  than  the  incident  ion 
energy.  Ion  implantation  is  done  at  ion  energies  where  the  sum  of  all  the  energy  lost  in  collisions  with 
electrons  and  the  sum  of  all  the  energy  lost  in  collisions  with  atoms  are  of  the  same  order-of- 
magnitude.  The  details  of  the  energy  transfer  of  ions  in  solids  is  a  large  subject  (Butler,  1981;  Borders, 
1980)  and  beyond  the  scope  of  this  chapter.  For  our  purposes  it  is  sufficient  to  note  that  the  Gaussian 
shape  of  the  concentration  profile  for  implanted  atoms  derives  from  the  statistical  distribution  of  possi¬ 
ble  energy  losses  in  random  and  independent  collisions  with  electrons  and  atoms  in  the  solid  as  the  ion 
slows  down.  Because  hundreds  to  thousands  of  collisions  typically  occur  before  the  ion  comes  to  rest, 
statistical  principles  may  be  applied  to  the  problem.  Indeed,  calculations  of  ion  range  parameters  are 
often  called  range  statistics. 

ii.  Ion  Range  Tables  for  Elemental  Targets 

A  number  of  tables  (Johnson,  1970;  Gibbons,  1985),  graphs  (Smith,  1977;  Winterbon,  1985), 
and  computer  codes  (Manning,  1974;  Biersack,  1980;  Biersack,  1983;  Manning,  1984)  may  be  found  in 
the  literature  which  predict  the  range  parameters  of  ions  in  solids.  We  have  included  in  this  section 
another  range  table  which  differs  from  those  previously  published  in  several  ways.  The  entries  in  Table 
(4.1)  were  calculated  using  a  modified  version  of  the  code  E-DEP-1  (Manning,  1974:  Manning,  1984). 
The  modifications  are,  i)  improved  estimates  for  the  atomic  collision  cross  section  (Wilson,  1977)  and, 
ii)  a  semi-empirical  adjustment  of  the  electronic  collision  cross-sections  (Land,  1978).  With  these 
modifications,  which  apply  to  all  ions  in  the  table  heavier  than  boron,  the  code  is  the  most  accurate 
available  for  a  Boltzmann  transport  approach  to  range  calculations.  Even  with  these  modifications  the 
tabulated  values  are  accurate  only  to  ±15%  in  Rp  and  ±20%  in  A Rp.  Relative  values  between  dif¬ 
ferent  ions  in  a  substrate  or  different  ion  energies  should  be  considerably  more  accurate.  The  tables  are 
simplified  and  more  concise  than  previous  tables  because  i)  ions  and  substrates  of  interest  for  metal 
implantations  are  emphasized,  ii)  only  two  or  three  significant  figures  are  included,  iii)  only  five  ion 
energies  (10,  50,  100,  200  and  400  keV)  per  ion/target  combination  are  included  because  linear  inter¬ 
polation  between  the  energies  is  accurate  to  a  few  percent  which  is  much  better  than  the  errors  in  the 
absolute  value,  and  iv)  the  values  of  Rp  and  A/?p  are  given  in  nm  rather  than  the  usual  mg/cm2 
Seventeen  ions  and  Seventeen  substrates  are  included  in  Table  (4.1).  The  usual  nomenclature  to 
define  the  atomic  numbers  and  atomic  weights  is  as  follows:  Zh  M\  are  the  atomic  number  and  the 
atomic  weight  of  the  ion,  respectively;  and  Z2  and  M2  are  the  atomic  number  and  atomic  weight  of  the 
target,  respectively. 

For  ions  other  than  those  in  the  Table  (4.1),  linear  interpolation  for  Rp  and  A Rp  using  the  atomic 
number  Z,  of  the  ions  is  accurate  to  better  than  4%.  For  targets  other  than  those  in  the  table,  one 
must  first  convert  the  appropriate  table  entries  to  m&'cm2  by  multiplying  by  the  density  of  the  targets, 
perform  the  linear  interpolation  with  the  mg/cm2  values  using  the  atomic  numbers,  and  than  reconvert 
the  answer  to  nm  by  dividing  by  the  density  of  the  elmental  target  of  interest. 


In  many  cases,  it  isn't  necessary  to  intepolate  (or  ions  or  targets.  If  the  ion  is  not  in  Table  (4.1) 
and  only  differs  from  the  nearest  ion  in  the  table  by  I  or  2  atomic  numbers,  than  the  values  of  the 
nearest  ion  will  be  adequate  Similarly,  for  targets  that  are  not  in  the  table  and  which  only  differ  from 
the  nearest  table  entry  by  I  or  2  atomic  units,  adequate  values  for  Rp  and  A/?;,  may  be  obtained  by 
scaling  with  the  mass  densities  of  the  unknown  target  and  the  table  entry.  For  example.  Rr  for  100 
keV  Cr  in  Ni  obtained  from  the  table  entries  for  Fe  is, 

Rp  (Cr  —\i )  -  Rp  (Cr  -  Fe)x^ f  =  38  x  =  34  nm 

f).\i  o.V 

which  agrees  with  the  table  entry  for  100  keV  Cr  into  Ni  For  target  atomic  numbers  Z:>  60  this 
method  can  be  used  with  good  accuracy  for  any  target  between  table  entries.  For  example,  for  100  keV' 
Zr  ions  we  have  Rp( Zi — -Ta)  =  Rp(Zr—Pi)  x  21.4/16.6  =  24nm  which  is  the  same  value  as  the  table 
entry  lor  Zr— * Ta.  A  list  of  the  densities  of  the  elements  is  given  in  Table  (4.2)  for  this  purpose  For 
steels,  the  table  values  for  iron  can  be  used,  scaled  by  the  density  of  steel  (e  g  .  (Rp  or  A Rp ) W(r/  = 
uri  x  < Rp  or \Rp  >fl  .),  The  densities  of  steels  range  from  7.5  to  8.9  g/cm '. 

Frequently,  the  target  is  an  alloy  composed  of  two  or  more  elements  There  are  two  methods  to 
derive  Rr  and  ±Rr  lor  these  cases  An  approximate  method  is  to  find  Rp  and  ±Rr  for  the  major  con¬ 
stituent  and  scale  the  results  by  the  ratio  of  the  density  of  the  element  in  the  table  and  the  density  of 
the  alloy  For  alloys  or  compounds  where  there  is  a  large  fraction  of  elements  of  largely  different 
masses  a  different  approximate  approach  is  required  The  table  entries  for  Rp  and  A Rp  must  be  con¬ 
verted  to  mg/ cm-  by  multiplying  by  the  respective  densities  and  a  weighted  sum  of  the  mg/cnv  values 
performed  The  result  is  reconverted  to  nm  by  dividing  by  the  density  of  the  alloy  or  compound  For 
an  alloy  or  compound  A.B,  where  x  +  y  =  l.  we  have. 


(4  1) 


(4.2) 


A  Rf"'  = 


A  R‘ 


(4  3 ) 


This  estimate  of  Rr  and  A Rr  for  alloys  is  an  approximation  because,  strictly  speaking,  the  calculation 
should  be  performed  by  including  different  stopping  powers  of  the  constituents  of  the  alloy  in  the 
Bolt/man  transport  equations  to  arrive  at  Rp  and  A Rp  However,  the  approximation  works  well  enough 
to  provide  a  1st  order  estimate  of  Rr  and  A Rp  For  determination  for  compounds  like  stoichometnc 
oxides,  carbides,  nitrides,  etc.  publications  mentioned  earlier  must  be  consulted  because  Table  <4  1) 
includes  entries  only  for  solid  elemental  targets 

Some  values  ol  A R.  are  not  listed  in  the  table  This  arises  because  the  numerical  calculation  of 
Johnson  and  (iibbons  range  parameters  on  which  the  code  F-DFP-I  is  based  is  unstable  lor  these  cases 
and  no  solution  exists 


Occupy  ing  onlv  two  book  pages,  this  is  the  most  compact  range  table  in  existence  and  allows  con¬ 
venient  estimates  of  Rf,  and  A Rr  lor  most  um/target  combinations  l  ittle  precision  is  lost  in  the  inter¬ 
polations  and  approximations  necessary  tor  arriving  at  values  for  ions  and  targets  not  included  in  the 
table  because  the  starting  values  ol  the  range  parameters  are  onlv  accurate  to  15"*. 


Table  4.2:  Selected  Properties  of  the  Elements 
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Number 

(j2> 


Symbol 


N 

(xlO22 
atoms  /cm3) 


(#/cm3) 


Atomic 

Number 

(r2) 


Symbol 


N 

(xlO22 

atoms/cm3) 


P 

(g/ cm) 


C.  Sputtering 


When  an  ion  strikes  and  penetrates  the  surface  of  a  solid,  some  of  the  near  surface  target  atoms 
receive  enough  kinetic  energy  in  the  collision  cascade  to  overcome  the  surface  binding  energy  and  are 
ejected  from  the  surface.  The  parameter  used  to  describe  this  phenomenon  is  the  sputtering  coeffi¬ 
cient,  S,  defined  as  the  ratio  of  the  number  of  atoms  ejected  from  the  surface  per  incident  ion.  For 
example,  values  of  S  for  100  keV  ions  incident  on  iron  vary  from  0.05  for  He  ions  to  13  for  Bi  ions. 
For  high  dose  implantation,  the  sputtering  coefficient  is  the  major  parameter  determining  the  depth  dis¬ 
tribution  and  surface  concentration  of  implanted  ions.  For  this  reason  sputtering  theory  is  presented  in 
enough  detail  to  enable  the  estimation  of  S  for  any  ion/target  combination. 

The  most  useful  theory  is  that  of  P.  Sigmund  (Sigmund,  1969)  who  proposed  the  following 
expression  for  the  sputtering  coefficient, 

s  =  0.042  a  S„  (eV/A>  1 

Ua(eV)  N (atoms/ $)  (cos0)5/} 

where  a  is  a  dimensionless  constant  and  is  a  function  of  M2/Mu  S„  is  the  atomic  stopping  power  (also 
called  nuclear  stopping)  computed  from  the  equations  which  follow  and  contains  the  dependence  of  S 
on  Zb  Z2,  M2  and  ion  energy  E.  U0  is  the  surface  binding  energy,  N  is  the  atomic  density,  and  0 
is  the  angle  of  incidence  to  the  surface  normal.  For  convenience,  values  of  U0  and  N  are  given  in 
Table  (4.2).  Tabulated  values  of  N  have  to  be  multiplied  by  10~24cm3/$  to  be  in  the  correct  units  for 
the  equation.  Entries  for  U0,  taken  from  the  Handbook  of  Chemistry  and  Physics  (CRC,  Cleveland 
1966),  are  the  "Heats  of  Formation  of  Gaseous  atoms  from  Elements  in  their  Standard  States"  con¬ 
verted  from  Kcal/mole  to  eV.  Because  of  approximations  used  in  the  theory,  equation  (4.4)  is  also  re¬ 
stricted  to  values  of  Mt>  14  where  the  electronic  stopping  power  is  a  small  fraction  of  the  nuclear 
stopping  power,  to  values  of  M2/Mx<4,  and  to  values  of  0<7O°  In  addition,  for  M\<  <M2 ,  S  a  (cos 
0)-1.  For  Mx<  14,  equation  (4.4)  overestimates  the  sputtering  coefficient.  For  ion  masses  greater 
than  A/,  =  120  incident  on  substrates  with  masses  greater  than  M2  =  120,  thermal  spike  effects  can 
sometimes  increase  S,  but  the  amount  of  the  increase  is  always  less  than  1.8  and  it  is  usually  much 
smaller. 

The  parameter  a  can  be  converted  to  an  approximate  analytical  form  for  convenient  calculation  of 
S  with  the  expression  (Matsunami  et  al.  1984), 


«  =  0.08  +  0.164 


M2 

~M~X 


0.4 


+  0.0145 


My 


M, 


1  29 


(4.5) 


Substituting  this  expression  into  equation  (4.4)  for  S,  and  noting  all  the  restrictions  we  have 
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Equation  (4  6)  indicates  that  S  is  directly  proportional  to  the  atomic  stopping  power  Sn  and  inversly 
proportional  to  the  sublimation  energy  U„.  The  angular  dependence  is  fairly  strong  indicating  that  S 
increases  by  factors  of  1.3  and  3.2  at  30°  and  60°,  respectively.  The  restriction  that  S  is  not  valid  for 
M |<  14  is  not  a  serious  problem  because  sputtering  coefficients  are  small  for  light  ions  and  the  range 
distribution  is  well  represented  by  the  Gaussian  distribution.  The  atomic  stopping  power  S„  is  com¬ 
puted  from  the  semi-empirical  expression  of  Wilson  et  al.  (Wilson,  1977)  because  it  is  more  accurate 
than  that  originally  used  by  Sigmund  and  because  of  convenient  formulation  The  atomic  stopping  is 
given  by. 
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Epsilon  is  the  reduced  energy  common  to  LSS  range  theory  (Lindhart,  1963)  but  with  the  Firsov 
screening  length  Za  instead  of  the  Thomas-Fermi  screening  length  as  prescribed  by  Wilson  et  al.  (Wil¬ 
son,  1977).  Sr  is  the  Wilson,  Haggmark  and  Biersack  atomic  stopping  power  in  reduced  energy  units 
and  equation  (4.7)  converts  S„  into  units  of  eV/A.  The  expression  for  S„  is  a  simple  algebraic  formula 
containing  Z,,  Z2,  M\,  M2 ,  the  ion  energy  E  and  mass  density  p  of  the  substrate  and  is  easily  pro¬ 
grammed  in  Fortran  or  Basic  for  the  purpose  of  computing  S.  The  values  of  S„  are  valid  for  all 
ion/target  combinations  listed  in  the  range  table  for  energies  between  10  and  400  keV. 

Equation  (4.6)  predicts  the  trends  of  sputtering  data  in  the  literature  quite  well  but  is  likely  to 
miss  the  absolute  value  of  S  by  as  much  as  a  factor  of  two.  Reynolds  and  Sigmund  (Reynolds  (1985),  ; 
Sigmund,  1969;  Sigmund,  1981)  have  discussed  the  reasons  for  these  discrepancies,  none  of  which  are 
easily  estimated.  The  two  largest  effects  are  ion  channeling  in  textured  polycrystalline  targets  and  sur¬ 
face  contamination  during  sputtering,  both  of  which  tend  to  reduce  S. 

Calculated  values  of  S  for  the  Group  IVB  and  VB  elements  such  as  Ti,  Nb  and  Ta  on  Fe  are  usu¬ 
ally  larger  than  experimental  measurements.  There  is  a  well  known  experimental  problem  in  measuring 
the  sputtering  coefficient  for  these  elements  because  of  their  reactive  nature.  Impurities  chemisorbed 
on  the  surface  can  greatly  influence  the  sublimation  energy  which  in  general  causes  the  experiment  to 
yield  a  value  smaller  than  is  calculated.  For  most  other  elements,  the  values  for  the  sputtering  coeffi¬ 
cient  S  should  predict  the  correct  value  to  within  25%.  The  best  choice  of  S  would  an  experimental 
value  from  the  literature.  Two  excellent  sources  of  values  of  S  are  references  Behrisch,  1981  and 
Matsunami,  1984.  If  S  is  known  for  a  given  system  but  for  a  different  energy,  then  the  experimental 
values  can  be  used  to  calibrate  equation  (4.6)  from  which  accurate  values  of  S  can  be  obtained  for  other 
ion  energies. 

D.  Atomic  Fraction  of  Implanted  Ions 

To  plan  a  particular  experiment  or  treatment  of  a  part,  it  is  important  to  be  able  to  predict  the 
concentration  profile  of  the  implanted  impurity  atoms.  This  section  discusses  the  factors  determining 
the  profile  in  low  dose  and  high  dose  regimes  and  gives  expressions,  tables  and  graphs  to  allow  predic¬ 
tions  of  the  impurity  profile  to  be  made. 

i  Low  Dose  Concentration  For  All  Ion  Masses 

For  low  ion  doses  at  normal  incidence  where  the  amount  of  material  sputtered  from  the  surface 
S6/ N  is  ^0  15/?,,,  the  ions  are  completely  retained  in  the  surface  and  the  concentration  profile  can  be 
approximated  by  a  Gaussian  distribution  characterized  by  a  centroid,  Rp,  a  standard  deviation.  ARP, 
and  a  peak  atomic  density,  Np,  given  as. 
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Assuming  the  target  density  doesn’t  change,  the  peak  atomic  fraction,  Cp ,  is. 


(4.11) 


(4.12) 


where  N  is  the  atomic  density  of  the  substrate  given  as  pA/ M2  where  A  is  Avogadro's  number  =  6.02 
x  1023  atoms/mole.  For  non-normal  angles  of  incidence,  Rp  is  shortened  by  a  cos  9  factor  so  that, 

Rp(9)  =  Rn  cos  9.  (4.13) 


Because  lateral  straggling  spreads  the  implanted  ion  distribution  in  a  direction  parallel  with  the  surface 
by  about  the  same  magnitude  as  the  straggling  A  Rp,  it  is  a  good  Ist-order  approximation  to  assume  that 
A Rp  is  independent  of  the  angle  of  incidence.  This  will  be  discussed  in  more  detail  in  the  following 
section. 


The  ion  dose  a  sample  receives  is  always  given  with  respect  to  an  area  parallel  to  the  sample  sur¬ 
face  but  the  measurement  of  dose  is  always  with  respect  to  an  area  normal  to  the  beam  direction. 
Therefore,  for  non-normal  angles  of  incidence  the  measured  dose  has  to  be  increased  by  a  cos  9  factor 
to  achieve  the  required  dose  in  the  area  projected  on  the  sample  surface,  or. 


meas. 
i  to  beam 


i  implanted 
>*11/0  surface 

cos  9 


(4.14) 


If  <t>  measured  stays  constant,  then  substitution  of  equation  (4.11)  and  (4.14)  into  equation  (4.12) 
shows  that  Cp(9)  =  Cp(9=  0)  cos  9.  This  is  the  angular  dependence  of  the  dose  in  the  low  dose  limit 
for  a  stationery  cylinder  when  the  ion  beam  floods  the  entire  diameter. 


Most  implantation  machines  have  a  minimum  energy  at  which  they  can  operate.  Implanting  at  a 
non-normal  angle  of  incidence  is  a  means  of  obtaining  shallow  implants  when  it  is  desired.  However, 
the  sputtering  coefficient  increases  sharply  with  angle  so  that  care  must  be  taken  that  the  condition  that 
the  amount  of  material  sputtered  away  be  <  0.15  Rp  still  applies. 

ii.  High  Dose  Concentration  for  Light  Ions:  M ]<?  20 


The  sputtering  coefficient  for  light  ions  (especially  H,  He,  Li,  B,C,  N  and  O)  is  generally  small 
enough  ( <  1 )  so  that  the  sputtering  can  be  ignored.  The  concentration  of  implanted  impurity  atoms 
builds  up  in  a  Gaussian  shape  at  depth  Rp  without  appreciable  erosion  of  the  surface.  In  this  case,  the 
atomic  fraction  at  the  peak  of  the  Gaussian  distribution  can  be  obtained  by  modifying  equation  (4.12) 
to  include  the  sum  of  the  target  density  N  and  the  implanted  atom  density,  or 


N„ 


Nn  +  N 


(4.15) 


iii  High  Dose  Concentration  for  Heavy  Ions:  A/, >20 

The  final  concentration  profile  in  the  high  dose  regime  is  the  result  of  many  processes  occurring 
simultaneously  These  processes  include  sputtering,  preferential  sputtering  of  the  surface  fraction  of 
one  element  when  the  target  is  an  alloy  or  a  compound,  radiation  enhanced  diffusion,  radiation  induced 
segregation,  ion  mixing,  compound  formation,  surface  reactions  with  gas  atoms  in  the  vacuum,  as  well 
as  other  processes  A  comprehensive  discussion  of  these  effects  is  beyond  the  scope  of  this  chapter 
(Borders.  1980)  Furthermore,  there  is  no  comprehensive  theory  that  can  predict  concentration  profiles 
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for  high  dose  implantations  because  of  the  interdependence  of  the  complicated  effects  listed  above.  It 
would  still  be  useful  however,  to  devise  a  means  to  give  a  Oth  order  approximation  to  the  profile  at 
high  doses  so  that  a  reasonable  estimate  of  ion  dose  can  be  made  prior  to  treatment.  To  that  end,  this 
section  describes  Oth  and  1st  order  approximation  methods  to  predict  the  saturation  dose,  the  surface 
concentration  of  implanted  atoms  and  the  concentration  profile  at  saturation  dose.  It  should  be  noted 
that  a  computer  code  is  under  development  that  includes  some  of  the  effects  listed  above  (Davisson. 
1986)  but  even  this  code  still  doesn't  include  all  1st  order  effects. 


For  high  ion  doses  (e  g.  >  101  ions/cm2)  removal  of  material  from  the  surface  by  sputtering 
uncovers  some  of  the  previously  implanted  ions  which  are  also  sputtered  from  the  surface.  As  a  result, 
less  than  1 00‘ ..  of  the  implanted  ions  are  retained.  This  effect  causes  large  deviations  from  the 
Gaussian-like  behavior  of  the  implanted  atom  profile.  A  Oth  order  approximation  for  the  surface  con¬ 
centration  at  saturation,  C\ur, ,  the  implantation  profile,  C(x),  and  the  saturation  dose.  <bs,  are.  respec¬ 
tively  (Liau.  1980) 
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where  .V,  is  the  atomic  fraction  of  implanted  atoms  and  S'B  is  the  atomic  fraction  of  substrate  atoms. 

Note  that  the  surface  concentration  is  the  same  as  the  bulk  concentration  and  the  bulk  concentra¬ 
tion  is  approximated  as  constant  from  the  surface  to  the  depth  Rn  Saturation  is  assumed  to  occur 
when  the  amount  of  material  sputtered  from  the  surface  reaches  the  depth  Rr .  e  g  .  S<£(  =  Rr\.  from 
which  equation  (4.18)  obtains.  Figure  (4.3)  illustrates  the  concentration  profile  obtained  for  these 
assumptions. 

A  different  approach  has  been  reported  by  Grabowski  et  al  and  Manning  (Grabowski.  1984;  Man¬ 
ning.  1985)  who  calculated  the  actual  concentration  profiles  under  certain  simplifying  assumptions. 
From  that  work  the  product  of  the  sputtei  coefficient  and  atomic  fraction  is  plotted  versus  a  dimension¬ 
less  depth  variable  x/(^/2\Rp)  for  the  cases  of  Rr  =  6  y/2^Rp  and  Rr  =  \/2±Rr  in  Figure  (4,4). 
These  curves  assume  the  implanted  ions  occupy  no  volume,  the  sputter  coefficient  and  Rr  remain  con¬ 
stant  during  implantation,  and  that  no  diffusion  processes  occur  The  latter  condition  could  be  realized 
it  the  implantation  were  performed  at  a  temperature  of  0°  K  These  assumptions  are  quite  restrictive 
but  provide  adequate  first  order  estimates  of  the  effects  of  sputtering  on  the  concentration  profile 
Sputter  saturation  is  said  to  occur  when  the  concentration  profile  no  longer  changes  with  increasing  ion 
dose  because  equilibrium  is  achieved  between  the  number  of  ions  being  implanted  and  the  number  of 
implanted  ions  being  removed  This  occurs  Corfu  =  3/2  in  Figure  (4  4a)  and  w  =  3  in  Figure  (4  4b) 
where  o<  is  the  amount  of  material  removed  relative  to  Rp  and  is  given  by 


In  Figure  *4  4a>  it  is  clear  that  for  <o  =  1/6  the  median  of  the  profile  is  only  slights  shifted  toward  the 
surface  by  sputtering  As  the  dose  is  increased  the  profile  fills  in  toward  the  surface  and  at  saturation 
the  concentration  profile  approaches  that  of  equation  (4  17)  In  figure  (4  4b).  which  represents  a  more 
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realistic  case  for  transition  elements  implanted  into  transition  metals,  there  is  no  well  defined  plateau  in  j 

the  profile  for  any  value  of  to.  The  value  of  to  for  the  Oth  order  estimate  for  saturation  dose  in  equa¬ 
tion  (4.19)  is  1.0,  whereas  the  saturation  values  of  to  in  Figures  (4.4a)  and  (4.4b)  are  1.5  and  3.0,  ) 

respectively.  This  shows  that  the  Oth  order  calculation  tends  to  underestimate  the  co,. 

From  a  family  of  curves  (Manning,  1985)  such  as  those  in  Figures  (4.4a)  and  (4.4b)  we  have 
developed  the  curves  in  Figure  (4.5)  for  to  at  the  saturation  dose  versus  Rp/ (V2A/?P)  to  allow  the  cal¬ 
culation  of  to ,  for  any  arbitrary  ion/ target  combination.  Because  the  surface  concentration  approaches 
the  saturation  value  assymptotically,  we  include  the  saturation  dose  which  achieves  90%  of  the  surface 
concentration  at  the  true  saturation  dose.  Otherwise,  considerably  more  ions  have  to  be  implanted  for 
the  limited  benefit  of  increasing  the  surface  concentration  by  10%.  A  curve  for  co  at  99%  saturation 
dose  if  also  shown.  Particularly  at  small  values  of  Rp/ (\fl\Rp)  the  dose  required  to  achieve  99%  j 

saturation  is  substantially  higher  than  for  90%  saturation.  Also  shown  in  Figure  (4.5)  is  the  parameter  I 

/3  which  is  a  correction  factor  applied  to  the  calculation  of  (-\url  ■  It  arises  from  the  fact  that  for  small 
values  of  Rpl  (\fl\Rp  ),  some  of  the  Gaussian  profile  resides  outside  the  surface  of  the  sample  and  is  . 

assumed  in  the  theory  to  be  lost  by  reflection.  Given  Rp ,  A Rp  and  S  one  can  find  $,  (90%)  by  solving 
equation  (4.19)  after  first  finding  to  (90%)  from  Figure  (4.5);  e.g.. 


(4.20) 


The  concentration  or  atomic  fraction  at  the  surface  for  the  $,(90%)  and  $,(99%)  are. 


$,(90%)  = 
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0.9/3 
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(4.21) 

Cur/  (99%)  = 

0.99/3 

S 

(4.22) 

where  it  is  assumed  that  there  is  no  change  in  the  substrate  density.  The  entire  profile  is  approximated 
by  equation  (4.17)  but  with  the  atomic  fraction  equal  to  C(x)  =  Cwr,  in  equation  (4.21)  or  (4  22). 


The  plot  of  (o  versus  Rp/\f2\Rp  in  Figure  (4.5)  may  be  accurately  represented  by  the  analytical 

form. 
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Substituting  equation  (4  23)  and  (4.24)  into  equation  (4.20)  then  gives  $,  for  normal  incidence; 

$ ,  ( 90% )  =  0.96  ~  (/?,  +  V2A/?,,)  (4  25) 

$,(99%)  =  ~  j/?p  +  1  64>/2AK„)  (4  26) 


lor  non-normal  angles  of  incidence  one  must  begin  with  the  corrected  values  of  Rp.  A Rp  and  S.  Sub¬ 
stituting  the  correct  angular  dependence  into  equation  (4  20)  we  have. 
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The  surface  atomic  fractipn  of  implanted  atoms  becomes. 


Qur/(90%)  =  (cos  9)511  ’  <4.28) 

These  expressions  show  that  as  the  angle  of  incidence  increases,  the  saturation  dose  and  surface  atomic 
fraction  decrease  by  the  factors  (cos  0)8/3  and  (cos  0)5/J,  respectively,  and  the  depth  of  the  implanted 
layer  decreases  by  the  factor  cos  0. 

The  angular  dependence  (cos0)8/3  for  <f>,  in  equation  (4.27)  is  a  good  1st  order  approximation  but 
it  is  not  explicity  correct  because  4>s  also  depends  upon  Inserting  equation  (4.23)  into  equation 
(4.27)  yields, 

<t>5  (90%)  =  0.96  ~  (cos  0)s/3  [/?,  cos  9  +  s/2A/?„]  .  (4.29) 

Note  that  only  if  A Rp(9)=bR„  cos  9  does  the  (cos  0)8/3  dependence  strictly  apply  which  may  be  shown 
by  factoring  cos  9  out  of  the  sum  in  equation  (4.29).  The  correct  angular  dependence  of  A Rp  is 

A/?„(0)  =  [a/?p2  (9  =  0)  cos2  9  +  A/?!2  (0  =  0)  sin2  0] '  2  (4.30) 

for  A/,  >  20  ;  0  <  45  0 

where  A/LL  is  the  lateral  straggling.  Lateral  straggling  refers  to  the  fact  that  ions  are  displaced  by 
atomic  collisions  from  their  initial  trajectory  so  that  the  fin.il  site  of  the  implanted  ions  has  a  statistical 
deviation  from  a  straight  line  path  described  by  A/LL,  or  straggling  in  a  direction  perpendicular  to  the 
initial  beam  direction.  Lateral  straggling  has  in  general  a  similar  magnitude  as  A Rp.  If  A/LL  =  A Rn 
then  the  straggling  is  isotropic  and  A  Rp  is  independent  of  the  angle  of  incidence.  We  have  computed 
A Rp(9)  from  equation  (4.30)  at  an  angle  of  45°  for  ions  ranging  from  B  to  Bi  and  substrates  from  Be  to 
Bi  in  order  to  assess  the  magnitude  of  the  angular  dependence  of  A Rp.  The  values  of  A/LL  were 
obtained  from  the  code  E-DEP-1  and  A  7^(0)  was  evaluated  for  energies  between  10  and  400  keV  that 
are  covered  in  Table  (4.1).  For  ion  masses  A/|  >  20,  A Rp(9)  deviated  from  A Rp  by  the  average  of 
-4%.  with  the  largest  deviations  being  -15%  and  +12%.  For  ions  with  A/,<20,  deviations  were  as 
high  as  +40%,  with  +20%  being  rather  typical.  Therefore,  equation  (4.29)  should  be  quite  accurate  if 
we  assume  that  A Rp  is  independent  of  angle  for  ion  masses  greater  than  about  20  and  for  angles  <  45°. 
The  deviations  of  A Rp{9)  from  A Rp  rapidly  become  less  as  the  angle  decreases  from  45°,  and  vice- 
versa.  For  ion  masses  less  than  20  the  sputtering  coefficient  is  frequently  small  enough  that  implanta¬ 
tion  to  saturation  is  not  practical  so  that  equation  (4.29)  wouldn’t  be  applicable  in  the  first  place. 

As  an  example,  consider  the  implantation  of  normally  incident  Cr  ions  into  Fe  at  an  energy  of  90 
keV.  Interpolating  Table  (4.1)  for  90  keV  we  have  Rp  =  38-0. 2x  (38-21)  =  35  nm  and  A Rp  =  18- 
0.2x(18-10)  =  16nm  which  compare  very  well  with  the  actual  E-DEP-1  values  of  34  and  16  nm, 
respectively.  The  sputtering  yield  has  been  measured  to  be  S  =  2.8  (S  differs  from  the  value  computed 
by  equation  (4.6)  because  of  texturing  of  the  polycryslalline  target)  (Reynolds,  1981,  1985).  From 
equation  (4.29)  we  have. 
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and  the  surface  concentration  from  equation  (4.21)  and  the  graph  in  Figure  (4.5)  for  (3  is. 


Cw,  (90%)  = 


0,9  x  0.99 
2.8 


0.32. 


For  an  angle  of  incidence  of  45°,  the  saturation  dose  computed  from  equation  <4  29)  becomes 
<*> ,  <90%)  =  y^-  (8  48  x  1022)  cos  (45  °js  *  [35  cos  45°  +  V2  16|  (lO  7| 

=  7.7  x  10lh  ions/cm2 


and  the  surface  concentration  becomes, 

C  (U„<90%)  =  °  9  *-°  89  (cos  45  V  1  =  0  16 

4.0 

where  f3  is  now  0.89  because  to  =  Rp/(J21Rr)  is  smaller  by  the  fact  that  the  range  ft,,  is  shortened  by 
a  cos  #  factor. 

iv.  Comparison  of  Calculated  Profiles  to  Experimental  Results 

Figure  (4.6a)  is  experimental  data  for  retained  dose  versus  angle  of  incidence  normalized  to  (he 
retained  dose  at  normal  incidence  for  several  ions  incident  on  steel  targets  and  for  doses  greater  than 
saturation  dose  (Grabowski.  1984).  The  retained  dose  should  also  to  first  order  follow  the  (cos  «)K  1 
dependence  The  agreement  is  good  for  angles  between  0  and  50°  which  indicates  that  the  simple  con¬ 
siderations  in  this  section  are  quite  acceptable  for  making  1st  order  estimates  of  <i,  and  C\url 
Discrepancies  can  arise  in  predicting  the  absolute  scale  of  the  (cos  #)"  ’  dependence  as  shown  in  Figure 
(4.6b)  for  Ta  implantation.  At  high  doses,  the  retained  dose  for  Ta-implanted  52100  steel  is  higher 
than  that  predicted  by  the  simple  theory  The  explanation  for  this  arises  from  two  effects  As  shown 
by  Grabowski  et  al.  (Grabowski,  1985),  diffusion  must  be  invoked  to  explain  the  observed  concentra¬ 
tion  profiles  of  Ta-implanted  iron  substrates.  Diffusion  of  Ta  into  the  bulk  would  tend  to  reduce  the 
surface  concentration  and  increase  the  retained  dose  compared  with  calculated  values  Secondly,  it  has 
already  been  mentioned  that  group  IVB  and  VB  elements  will  form  surface  carbides  during  implantation 
in  high  vacuum  (>10  Torr)  to  such  an  extent  that  atomic  concentrations  of  C'  up  to  15  at  %  are  not 
uncommon  (Singer,  1984)  Our  simple  theory  assumes  that  the  sputtering  coefficient  doesn't  change 
during  implantation  However,  carbides  in  general  have  a  lower  sputtering  coefficient  then  ihe  free 
metal  so  that  S  in  equation  (4.25)  decreases  as  the  carbon  content  of  the  surface  increases  This  el  feet 
tends  to  increase  the  retained  dose  and  the  surface  fraction  of  the  implanted  element  Figure  (4  7) 
demonstrates  this  effect  for  Ta-implanted  52100  steel  (Hubler,  1985).  Sputter  Auger  profiles  of  C  and 
Ta  for  different  doses  clearly  show  that  carbon  enters  the  surface  during  implantation  Therefore,  the 
calculation  of  <*>,  and  C\url  by  means  of  equations  (4,25)  and  (4  21)  should  be  considered  lower  limits 
for  the  refractory  carbide  forming  elements.  Singer  (Singer,  1984)  and  Behrisch  (Behrisch.  1981)  dis¬ 
cuss  this  problem  and  point  out  that  many  measurements  of  S  for  these  elements  in  the  literature  are 
effected  by  this  phenomenon  and  as  a  result  are  factor  of  2  to  3  too  low. 

A  definitive  experiment  demonstrating  the  reduction  in  sputtering  coefficient  with  surface 
absorbed  carbon  has  been  performed  by  Sartwell  and  Baldwin  (1985).  Ti  was  implanted  into  thin  bilayer 
films  of  iron  on  silicon  in  ultrahigh  vacuum  GO'9  Torr)  or  with  the  chamber  backfilled  to  a  partial 
pressure  of  10°  torr  with  CO.  The  total  amount  of  Ti  retained  in  the  sample  and  the  amount  of  Fe 
removed  from  the  film  were  continuously  monitored  by  means  of  proton  induced  X-ray  y  ields  Figure 
(4.8  a&b)  shows  the  X-ray  yields  vs  dose.  It  is  clear  that  in  the  backfill  case  the  sputtering  coefficient 
decreases  with  respect  to  the  UHV  implantation  and  that  more  Ti  is  retained.  The  sputtering  coeffi¬ 
cient  starts  at  a  value  of  3.5  finishes  at  a  value  of  1.8.  The  value  calculated  from  equation  (4  6)  is  4.1, 
in  reasonable  agreement  with  the  UHV  measured  value.  Therefore,  the  calculated  value  represents  an 
upper  limit  for  reactive  implanted  ions. 

Comparisons  of  calculated  and  experimental  values  measured  by  Reynolds  et  al  (Rey  nolds.  1981 ) 
of  <*>,  and  C\url  are  given  in  Table  (4.3).  The  experimental  points  were  generated  by  monitoring  the 
optical  signal  from  Cr  on  the  surface  of  iron  as  a  function  of  dose  during  Cr  implantation  The  optical 
signal  from  Cr  was  initially  zero,  increased  rapidly,  and  then  saturated  as  the  dose  increased  indicating 
that  <*> ,  had  been  reached  Similar  data  was  produced  for  Cu  implanted  into  Al  The  agreement  for 
C\url  and  <*>,(90%)  between  experiment  and  the  1st  order  values  is  quite  good. 


Table  4.3  —  Experimental  and  Calculated  Values  of  Csur)  and  <f>s  (90%) 


Measurement 
(Reynolds,  1981) 

0th  Order 

(equations  4.16,  4.18) 

1st  Order 

(equations  4.21,  4.25) 

r 

s urt 

0.24 

0.28 

0.25 

Cu  —  A1 

<t> ,  (90%) 

1.8  x  10!W 

1.1  x  1017/cm2 

1.6  x  10 l7/cm2 

90  keV 

(~\urt 

0.33 

0.36 

0.32 

Cr  —  Fe 

<f>,  (90%) 

1.7  x  10 '7cm2 

1.0  x  10n/cm2 

1.6  x  1017/ cm2 

90  keV 

E.  Selected  Experimental  Range  Profiles 

Some  of  the  most  widely  used  ions  for  the  implantation  of  steels  are  N,  Ti  and  Cr  so  considerable 
effort  has  been  expended  to  determine  concentration  profiles.  Examples  of  depth  profiles  for  these  and 
other  cases  are  given  here  in  Figures  (4.9)  through  (4.13).  In  Figure  (4.9)  the  nitrogen  concentration 
profile  departs  from  Gaussian-like  behavior  for  doses  greater  than  2  x  10 n/cm2  at  40  keV.  In  this  case, 
it  is  unclear  as  to  whether  radiation  enhanced  diffusion  or  sputtering  causes  the  increase  in  surface  con¬ 
centration  of  nitrogen  and  it  is  probably  caused  partly  by  both.  For  Nb  implanted  into  Fe  (Figure  4.12) 
there  is  no  difference  in  the  profile  for  doses  of  2.5  and  4.0  x  1017/cm2  showing  that  the  equilibrium 
between  implanted  ions  and  sputtered  ions  was  achieved  near  the  dose  of  2.5  x  10 nion/cm2.  The 
ranges  obtained  from  Table  (4.1)  are  in  reasonable  agreement  with  the  low  dose  data  in  Figures  4.7, 
4  9-4.12.  The  range  for  100  keV  Be  in  Ni  from  Table  (4.1)  is  1440  A  which  is  ~35%  smaller  that  the 
data  in  Figures  (4.13)  indicates.  This  large  discrepancy  arises  because  the  semi-empirical  corrections  to 
electronic  stopping  power  used  in  the  calculation  of  ranges  in  the  program  E-DEP-1  are  not  available 
for  elements  lighter  than  carbon.  An  overestimate  of  the  electronic  stopping  power  causes  the  com¬ 
puted  range  to  be  too  small.  The  independence  of  the  shape  of  the  concentration  profile  with  ion  dose 
for  light  ions  that  results  from  a  small  value  of  S  is  also  demonstrated  by  the  data  for  100  keV  Be  in 
Figure  (4.13).  Note  that  the  peak  position  of  the  profile  doesn't  change  for  doses  of  1016  and 
10  r  ions/ cm2. 

F.  Influence  of  Masking 

* 

A  common  problem  to  be  addressed  in  processing  is  the  implantation  of  cylindrical  or  spherical 
objects  By  masking  the  samples  to  allow  the  beam  to  impinge  only  at  angles  near  normal  incidence, 
concentration  profiles  can  be  calculated  by  the  simple  expression  used  above.  Opening  the  mask  to 
larger  angles  or  completely  flooding  the  sample  makes  more  efficient  use  of  the  ion  beam.  In  this  case 
the  depth  profiles  must  be  integrated  over  all  angles  of  incidence.  Figure  (4.14)  is  such  a  calculation 
for  the  case  of  B  implanted  into  Be  (Smidt,  1985).  Profiles  for  the  same  dose  are  shown  for  four  dif¬ 
ferent  masking  angles.  The  sputtering  coefficient  is  very  small  for  this  case  so  most  of  the  differences 
are  caused  by  range  shortening  with  angle  only.  More  dramatic  effects  that  can  occur  for  larger  sputter¬ 
ing  coefficients  are  shown  by  the  data  of  Singer  and  Jeffries  (Singer,  1984)  in  Figure  (4.15).  Profiles 
of  Ti-implanted  steel  are  measured  for  saturation  doses  implanted  into  52100  steel  with  3  different 
experimental  configurations. 


G.  Topographical  Changes 


m 


In  general,  the  larger  the  sputtering  coefficient,  the  greater  the  possibility  that  the  surface  will  be 
roughened  on  a  microscopic  scale  during  implantation.  Examples  of  roughened  surfaces  are  presented 
in  Figure  (4.16)  which  are  surface  profiles  of  AIS1  M50,  CRB7  and  440C  steel  3/8"  diameter  rods  after 
implantation  with  Ta  ions  (Hubler).  The  rods  were  rotated  during  treatment  and  a  mask  restricted  the 
angles  of  incidence  to  ±  30°.  The  region  profiled  spans  unimplanted  and  implanted  material  and  indi¬ 
cates  that,  i)  1800  Aof  material  was  removed  by  sputtering  for  M50  and  CRB7  and  1200  Afor  440C,  ii) 
the  surface  in  the  implanted  region  has  been  roughened  with  respect  to  unimplanted  material,  especially 
for  CRB7  and  440C.  Figure  (4.17)  presents  differential  interference  contrast  photomicrographs  of  the 
same  rods,  additional  unimplanted  rods,  and  rods  implanted  with  N  ions  in  the  same  geometry.  N 
implantation  has  little  effect  on  surface  finish.  Mesas  have  developed  on  the  Ta-implanted  rods  which 
result  from  a  faster  erosion  rate  for  the  iron  matrix  than  for  the  carbides  which  form  the  mesas.  The 
mesas  for  CRB7  appear  as  depressions  because  of  a  phase  reversal  of  the  interference.  This  roughening 
effect  is  minimized  by  i)  confining  the  beam  to  near-normal  angles  of  incidence  in  order  to  minimize 
the  sputtering  coefficient,  ii)  selecting  the  minimum  ion  dose  necessary  to  treat  the  sample,  and  iii)  by 
purposely  backfilling  the  target  chamber  with  a  hydrocarbon  (e.g.  CO,  C//4)  in  order  to  produce  surface 
carbide  which  reduces  the  sputter  coefficient  of  the  matrix  (Sartwell  and  Baldwin,  1985).  Only  strong 
carbide-forming  substrates  and  ions  are  amenable  to  the  latter  treatment.  Surface  roughening  is  very 
material  dependent  as  shown  by  comparing  the  micrographs  in  Fig.  (4.16)  for  M50  and  CRB7 
implanted  with  Ta  ions. 

Figure  (4.18)  compares  the  surface  topography  of  two  440C  steel  rods  implanted  with  Ta  ions  to 
identical  doses,  but  for  one  the  vacuum  was  1  x  10-6  Torr  and  for  the  other  the  target  chamber  was 
backfilled  to  a  pressure  of  5  x  10“'  torr  with  CO  gas.  Note  that  the  amount  of  material  removed 
decreased  by  a  factor  of  almost  three  and  the  surface  retained  its  initial  smoothness  in  CO  backfilled 


Surfaces  can  also  be  smoothed  by  ion  implantation.  For  example,  Pronko  et  al.  (Pronko,  1984) 
have  reported  surface  smoothing  of  polished  Mo  surfaces  by  Mo  ion  implantation.  Also,  removal  of 
very  fine  final  polishing  scratches  on  M50  and  52100  steel  is  a  frequently  observed  phenomenon. 

H.  Workpiece  Manipulation 

Tools  and  components  of  mechanical  machinery  can  have  widely  varying  shapes  and  complex 
geometries  which  present  specific  problems  for  ion  implantation  treatment.  Furthermore,  beam  pro¬ 
cessing  with  laser,  electron  or  ion  beams  is  line-of-sight  processing  which  places  some  limitations  on  the 
types  of  components  that  can  be  treated.  For  example,  the  internal  surfaces  of  small  tubes  or  holes 
cannot  be  ion  implanted. 

The  components  to  be  implanted  are  described  as  having  planar,  cylindrical  or  spherical  surfaces. 
Planar  surfaces  are  most  easily  implanted  because  they  frequently  require  no  manipulation  in  front  of 
the  ion  beam  or  only  some  lateral  translation  in  order  to  treat  the  entire  area.  Cylindrical  surfaces  are 
implanted  by  rotating  the  piece  in  front  of  the  ion  beam.  Long  cylindrical  objects  may  be  treated  by 
translating  the  rotating  cylinder  along  the  direction  of  the  cylinder  axis.  Figure  (4.19)  schematically 
shows  these  geometries,  including  the  geometry  for  treatment  of  the  inner  surface  of  a  ring  such  as  a 
bearing  race. 

I.  Workpiece  Cooling 

One  of  the  greatest  difficulties  in  high  dose  ion  implantation  processing  is  the  maintenance  of  low 
workpiece  temperatures.  For  a  beam  current  of  3  mA  at  an  energy  of  100  keV  the  incident  power  is 
300  watts.  This  power  must  be  dissipated  or  the  workpiece  will  be  overheated,  especially  for  carbon 
steels  whose  tempering  temperatures  may  be  as  low  as  125°C.  To  address  this  important  problem  we 
refer  to  the  work  of  Grabowski  and  Kant  (Grabowski.  1983)  who  have  reported  methods  to  remove 
heat  in  a  vacuum.  The  following  paragraphs  summarize  their  findings. 


In  a  vacuum,  heat  can  be  removed  by  radiation  or  by  conduction  through  solids  to  a  heat  sink. 
For  low  doses  or  large  mass  objects,  the  heat  capacity  of  the  object  can  be  sufficient  to  maintain  the 
temperature  under  a  specified  limit.  Radiative  cooling  is  effective  only  for  high  target  temperatures 
(>400°C)  or  large  surface  areas,  both  of  which  are  usually  not  realizable  for  high  dose  ion  implanta¬ 
tion  of  metals.  Figure  (4.20)  shows  the  power  radiated  per  unit  surface  area  as  a  function  of  surface 
temperature  for  emissivities  of  1.0  and  0.2,  the  latter  being  a  typical  value  for  steels.  The  dashed  lines 
show  that  at  temperatures  of  500°C  and  20 0°C  the  radiated  power  is  0.4  W/cm 1  and  0.05  W/cm 2, 
respectively.  Thus,  for  an  incident  power  of  300  Watts  a  surface  area  of  750  cm 2  is  required  to  main¬ 
tain  the  temperature  at  500°C  by  radiation  alone.  This  is  an  unacceptable  temperature  rise  for  many 
metals  and  the  surface  area  required  to  maintain  500°C  implies  a  part  which  is  near  the  maximum  size 
that  can  be  conveniently  placed  in  implantation  systems.  Moreover,  the  surface  area  must  increase  by  a 
factor  of  8  in  order  to  maintain  the  temperature  rise  to  less  than  200°C.  Another  means  to  use  radia¬ 
tive  cooling  is  to  reduce  the  beam  power  which  can  lead  to  unacceptably  long  implantation  times. 


Alternative  ways  to  make  use  of  high  beam  currents  are  by  conduction  through  the  workpiece  to  a 
heat-sink,  or  by  use  of  the  heat  capacity  of  the  part  as  the  heat-sink.  Grabowski  and  Kant  (Grabowski, 
1983)  have  analyzed  the  situations  of  one-dimensional  heat  flow  through  a  flat  plate,  a  cylinder  or  a 
sphere  and  the  results  are  plotted  in  Figure  (4.21).  Radiative  cooling  is  not  included  in  the  curves  so 
they  do  not  apply  for  temperatures  >500°C.  A  dimensionless  surface  temperature  rise 
Bs  =(  Ts—  Tc, )/ ( Jl/k )  is  plotted  versus  dimensionless  time  r=a<j>E/ l2J .  In  these  expressions  Ts  is  the 
surface  temperature,  T0  the  initial  temperature  and  Tb  is  the  temperature  either  at  the  back  of  a  plate 
of  thickness  I  or  at  the  center  of  a  sphere  or  cylinder  of  radius  1.  Tb  will  be  discussed  later  in  relation 
to  the  three  regimes  shown  in  the  figure.  The  other  variables  are  the  incident  power  density 
J {W/cm2),  the  thermal  conductivity  k(W/(cm-K))  the  exposure  time  t(s),  and  the  thermal  diffusivity 
aicm2/ s)  =  k/pC  where  p  is  density  ( g/cm3 )  and  C  is  heat  capacity  (Ws/(g-K).  </>  is  the  ion  dose 
(ions/ cm2)  and  E  is  the  ion  Energy  (Joules/ion).  Representative  values  of  a,  k  and  C  are  given  in 
Table  (4.4)  for  several  materials  of  interest. 


Table  4.4  —  Thermal  Properties  of  Selected  Materials  (from  Grabowski  and  Kant,  1983) 


Material 

(ctrr/s) 

k 

(W/(cmK)) 

PC 

(Ws/(cm 

Cu 

1.1 

3.9 

3.4 

A1 

0.95 

2.3 

2.4 

Al  (2024-T4) 

0.50 

1.2 

2.4 

Si 

0.53 

0.84 

1.6 

Brass  (Cu  -30Zn) 

0.38 

1.2 

3.2 

Ta 

0.24 

0.54 

2.3 

Fe 

0.23 

0.80 

3.5 

Fe  (0.100 

0.18 

0.65 

3.5 

Fe  (M2  tool  steel) 

0.058 

0.21 

3.7 

Fe  (304  SS) 

0.041 

0.16 

4.0 

Ti 

0.068 

0.16 

2.3 

Ti  (6A1  -  4V) 

0.026 

0.068 

2.6 

Superalloy  (IN  738) 

0.035 

0.12 

3.4 

Three  different  regimes  of  conduction  heat  transfer  are  indicated  in  Figure  (4.2).  In  regions  1  and 
II  all  deposited  energy  is  retained  in  the  sample.  In  region  III  a  steady-state  condition  is  reached  where 
heat  flow  is  constant  through  a  plate  sample  to  a  heat  sink.  In  region  I,  Tb  =  T„  because  the  deposited 
energy  has  not  yet  diffused  to  the  boundary  of  the  sample.  For  cylinders  or  spheres,  Th  is  at  the  center 
of  the  sample.  In  region  II  most  of  the  deposited  energy  is  distributed  throughout  the  sample  so  that 
Th^=Ts.  The  transition  from  region  I  to  either  II  or  III  occurs  at  r  a  1.  For  most  metal  implantations 
with  high  current  beams,  r  is  greater  than  1.  For  plate  samples,  expressions  for  the  temperature  rise  in 
regions  I,  II,  and  III  are. 
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11  T~T0 


(4.32) 


n  <(>  E 
P  C  I 


n  —  1  plate 
n  =  2  cylinder 
n  =  3  sphere 


111  Ts-T0  =  —  (ideal  heat  sink)  . 


(4.33) 


In  Figure  (4.22)  the  generalized  curves  of  Figure  (4.21)  are  evaluated  for  the  special  cases  of  a 
tool  steel  and  copper  plates  for  an  input  power  density  of  100  W/cm 2  (100  keV,  1  mA/cm 2)  and  a 
thickness  of  1  cm.  These  represent  two  extreme  cases  for  implantation  of  metals  because  steel  is  a 
relatively  poor  conductor  of  heat  whereas  copper  is  a  good  conductor.  In  10  seconds,  the  copper  and 
steel  surface  reach  temperatures  of  290°C  and  430°C,  respectively.  For  a  1  mA/cm 2  beam  it  requires 
16  seconds  to  implant  a  dose  of  1017/cm2.  Therefore,  the  steel  must  be  actively  cooled  in  order  to 
maintain  the  surface  temperature  below  the  typical  tempering  temperatures  for  carbon  steels.  Even  for 
the  case  of  conduction  to  an  ideal  heat  sink,  the  steel  reaches  a  steady-state  value  of  450°C.  Therefore, 
the  beam  power  must  be  reduced  or  a  duty  cycle  introduced  which  lowers  the  average  incident  power. 

Equation  (4.33)  for  region  III  implies  perfect  conduction  of  heat  from  the  sample  to  the  heat-sink 
which  in  practice  is  rarely  achieved  in  vacuum  because  temperature  drops  occur  across  interfaces  where 
the  sample  is  clamped  to  the  heat-sink.  The  temperature  differential  across  the  interface  is  expressed 
as, 

A  T-J/h  (4.34) 

where  h  ( W/cm2-K)  is  the  interface  contact  conductance.  According  to  Grabowski  and  Kant  (Gra- 
bowski,  1983),  h  can  be  maximized  by  using  a  pressure  contact  and  by  inserting  a  thin  foil  of  soft  high 
thermal  conductivity  material  at  the  interface  (e.g.  Cu,  In,  Al).  Table  (4.5)  presents  values  of  the  heat 
conductance  across  interfaces  in  vacuum  for  different  metal-metal  interface  configurations  measured  for 
a  contact  pressure  of  145  psi.  Table  (4.5)  provides  useful  guidelines  for  maximizing  the  interface  con¬ 
tact  conductance  for  most  situations  that  are  encountered  for  implantation  of  metals. 

Hirvonen  (1986)  has  analyzed  the  situation  where  a  cylinder  composed  of  tool  steel  is  cooled  by 
gripping  it  at  one  end  over  a  length  equal  to  1  diameter,  for  an  incident  power  density  of  1  W/cm2,  for 
a  thermal  conductivity  k=0.2  W/(cm-K),  and  for  an  interface  conductance  at  the  heat  sunk  end  equal 
to  1  W/(cm2- K).  Figure  (4.23)  is  as  plot  of  the  maximum  temperature  rise  of  the  tip  of  the  tool 
versus  the  tool  diameter  and  tool  length.  By  careful  interpolation  one  can  predict  the  maximum  tem¬ 
perature  rise  for  tools  ranging  from  0.01  inches  to  10  inches  in  diameter  for  lengths  between  1  and  10 
inches.  Examples  of  the  use  of  equation  (4.31-4.34)  are  given  in  the  final  section  V. 

J.  Examples  of  Fixtures  for  Manipulation  and  Cooling 

An  example  of  a  rotating  water-cooled  feed  through  used  to  cool  and  rotate  cylindrical  geometry 
samples  is  shown  in  Figure  (4.24)  (Grabowski,  1985a).  Proceeding  from  the  top  of  the  feed-through 
assembly,  there  is  a  rotating-seal  water  jacket  with  input  and  output  water  lines,  a  gear  to  attach  the 
motor  drive,  and  a  collar  around  the  4-cm-diameter  stainless  steel  shaft  to  hold  the  vertical  position  of 
the  shaft  against  the  actual  vacuum  feed-through.  The  vacuum  feed-through  is  sealed  by  a  ferrofluid  to 
provide  vacuum-tight  shaft  rotation.  It  is  mounted  on  a  welded-type  stainless  steel  bellows  with  4" 
diameter  flanges  to  allow  for  up  to  a  30  degree  tilt  angle  of  the  rotation  axis.  The  bellows  and  the 
mounting  assembly  which  secures  the  entire  feed-through  are  insulated  so  that  charge  may  be  collected 
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Table  4.5  —  Interface  Contact  Conductance 
(from  Grabowski  and  Kant,  1983) 


Interface 


hi W/(cm  -K)) 


ss/ss 

0.04 

SS/AI 

0.19 

SS/AI  Foil/SS 

>0.3 

SS/Indium/SS 

>0  6 

SS/Silver  Paint/SS 

>5 

Al/Brass 

0.5 

Al/Cu 

0.9 

AI/OF-CU 

1.5 

Al/AI 

16 

Brass/ Brass 

0  1 

Cu/Cu 

0.6 

OF-Cu/OF-Cu 

1.5 

SS:  304  Stainles  Steel 

Cu:  Tough-pitch  Cu 

Al:  6061 ,  T6  Temper 

OF-Cu:  Oxygen-free  Cu 

Brass:  Cu-30  Zn 

from  the  target  Jigs  which  hold  the  workpiece  are  secured  to  a  copper  plug  in  the  end  of  the  water 
cooled  feed-through  shaft  by  four  screws.  A  small  diameter  tube  inside  the  65-cm-long  hollow  shaft 
directs  water  at  the  copper  plug,  so  that  heat  generated  by  the  ion  beam  conducts  through  the  jig  and 
into  the  cooled  copper  plug  Dual  O-ring  seals  in  the  inner  diameter  of  the  ferrofluid  vacuum  feed¬ 
through  allow  linear  translation  of  the  shaft  without  breaking  the  vacuum  seal. 

Flat  plate  geometry  samples  can  be  secured  to  a  water-cooled  block  of  copper  or  aluminum  Fre¬ 
quently  used  methods  to  ensure  heat  conduction  from  sample  to  plate  include  pressure  from  a  clamping 
system,  silver  print  (organic  carrier  with  Ag  dust  in  suspension),  commercially  available  heat-sink  com¬ 
pounds  li  e.  MUNG  II,  Commonwealth  Scientific),  and  high  vacuum  fluids  (i.e.  Santovac  5)  with  Ag 
powder  in  suspension.  For  batch  processing  of  parts  a  duty  cycle  can  be  introduced  by  mechanical 
manipulation  in  order  to  lower  the  average  beam  power  on  individual  parts.  This  places  loss  stringent 
cooling  requirements  on  the  fixture  that  manipulates  parts  without  sacrificing  through-put  (Smidt  and 
Sartwell.  1985)  In  high  current  ion  implanters  (beam  currents  >  I  mA>  space  charge  blow-up  of  the 
ion  beam  precludes  the  use  of  a  raster  scanned,  focussed  ion  beam  so  mechanical  manipulation  or 
defocussing  are  required.  One  way  to  address  the  problem  is  to  develop  a  defocussed  ion  beam 
(Armini  and  Bunker,  1985)  which  is  capable  of  projecting  a  10  mA  beam  into  a  40  cm  diameter  area 
with  good  uniformity  of  beam  density  within  the  diameter.  In  this  scheme,  the  average  beam  current 
density  on  the  samples  is  reduced  to  10 fxA/cm},  a  value  at  which  the  removal  of  heat  from  the  work- 
piece  is  more  easily  managed  Figure  <4.25)  is  a  photograph  of  this  apparatus  and  Figure  <4  26)  shows  a 
commercially  available  non-mass  analyzed  ion-implanter  (Hirvonen;  1984)  that  uses  a  low  melting  point 
(145°C  )  eutectic  metal  ( 60%  pb/40"/<.  Bi)  to  both  hold  tools  and  act  as  a  heat-sink  medium. 

K.  Cleaning  of  Workpiece 

As  with  all  coating  techniques,  surface  cleanliness  prior  to  implantation  is  essential.  In  the  case  of 
coatings,  cleanliness  promotes  good  adhesion  of  the  film  to  the  substrate.  There  is  no  adhesion  prob¬ 
lem  with  ion  implantation,  but  cleanliness  ensures  a  uniform  implantation  dose  on  all  areas  scanned  by 
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the  beam  Since  the  ions  do  not  penetrate  dust  particles  or  oil  films  that  are  more  than  several  hundred 
nanometers  in  thickness,  the  areas  beneath  will  not  be  implanted  Cleanliness  and  surface  finish  are 
most  important  for  corrosion  and  oxidation  applications  where  unimplanted  regions  may  corrode  at 
rates  several  orders-of-magnitude  faster  than  implanted  material.  These  weak  areas  act  as  initiation 
sites  where  corrosion  will  penetrate  deeply  or  undermine  the  implanted  area.  If  the  surface  is  pitted 
prior  to  implantation  such  that  internal  surfaces  are  present,  the  ion  beam  cannot  reach  these  areas  and 
corrosion  may  proceed  unchecked  in  the  pits 

Both  cleanliness  and  surface  pits  are  not  as  important  for  wear  applications  where  occasional  small 
unimplanted  regions  represent  only  a  small  fraction  of  the  total  area  to  which  the  load  is  applied  How¬ 
ever.  it  is  always  generally  true  for  ion  implantation  or  thin  film  coatings  that  "cleanliness  is  next  to 
Godliness"  The  primary  culprits  to  be  aware  of  are  oil-and  grease-based  preservative  treatments  used 
on  many  steel  parts  to  keep  them  from  rusting,  oil  films  on  cold-rolled  sheet,  fingerprints  and  dust  par¬ 
ticles  It  is  good  practice  to  clean  metal  parts  with  solvents  whether  or  not  oil  films  are  known  to  be 
present  Methods  frequently  used  are  soaking  in  baths  of  paint  thinner,  toluene  and  high  purity 
methanol  or  ethanol,  in  that  order  Ultrasonic  action  improves  the  cleaning  when  it  is  compatible  with 
the  part  being  cleaned  (e  g  .  do  not  ultrasonically  clean  A)  or  Cu  because  they  damage  by  cavitation. 
Degreasing  in  freon  vapor  is  frequently  substituted  for  the  paint  thinner  step. 

The  best  practice  is  to  clean  the  workpiece  immediately  prior  to  performing  the  implantation  and 
to  blow  a  jet  of  oil-free  air  or  dry  nitrogen  over  the  surface  just  before  insertion  into  the  vacuum 
chamber  to  remove  dust  particles  If  this  procedure  cannot  be  used  then  storage  of  the  cleaned  work- 
piece  in  a  dessicating  container  for  long  periods  or  a  laminar  flow  bench  for  short  periods  is  recom¬ 
mended  Plastic  disposable  gloves  should  be  used  at  all  times  to  avoid  finger  prints.  In  short,  any 
available  means  of  providing  a  clean  environment  throughout  the  handling  of  the  workpiece  should  be 
employed  and  care  should  be  taken  to  ensure  that  the  surfaces  to  be  implanted  are  free  from  contami¬ 
nation  and  flaws  such  as  pits. 

L.  Size  Of  W  orkpiece 

The  si/e  of  a  workpiece  that  may  be  implanted  depends  upon  the  kind  of  vacuum  end-station  that 
is  peculiar  to  each  vendor  and  whether  the  ion  required  is  a  solid  or  a  gas  in  its  elemental  form.  For 
ions  capable  of  being  formed  from  a  gaseous  element,  a  non-mass  analyzed  beam  may  be  used  which 
lends  itself  to  the  implantation  of  larger  surface  areas  than  do  mass  analyzed  ion  beams  However,  new 
implantation  machines  with  mass  analysis  capable  of  10  mA  beam  current  are  becoming  available  so 
that  this  distinction  is  likely  to  disappear  in  the  near  future 

Because  of  differences  in  the  size  of  the  implantation  vacuum  chambers  employed  by  different 
vendors,  a  useful  way  to  describe  the  size  of  workpieces  is  to  give  examples  of  items  that  have  been 
implanted  A  partial  list  includes:  4-cylinder  automobile  engine  crankshafts  (Sioshansi.  1480).  10" 
diameter  H)  D  )  bearing  races  and  7/8"  diameter  balls  (48  balls/batch)  (Grabowski.  1985.  Popgoshev. 
I98.li.  thermo-plastic  injection  molds  and  feed  screws  (Dearnaley,  1983),  steel  pressing  tools  and  WC' 
punches,  dies  and  drills  (Sioshansi.  1985.  Dearnaley,  1983;  Hirvonen.  1984)  Batch  processing  of  a 
large  number  of  small  pieces  such  as  drills  is  straightforward.  Approximate  maximum  areas  for  implan¬ 
tation  are  12"  *  12”  for  planar  geometries  and  6"  x  30"  for  cylindrical  geometries  where  6"  is  the  diam¬ 
eter  The  group  at  the  Atomic  Energy  Research  Establishment  at  Harwell,  U  K.  has  a  non-analyzed- 
beam  ion  implanter  which  can  treat  an  area  I  meter  by  i  meter  (Dearnaley,  1984). 

VI.  Implantation  Time  and  Cost 

The  time  necessary  to  implant  a  sample  depends  upon  the  beam  current  available  at  the  vendor's 
facility,  the  maximum  allowable  temperature  for  the  material,  and  the  sample  geometry.  The  latter 
determines  the  beam  utilization  efficiency  Smidt  and  Sartwell  (Smidt.  1985)  have  analyzed  the  costs 
associated  with  running  a  Manufacturing  Technology  Facility.  Table  (4.6)  is  extracted  from  their  work 
and  shows  the  projected  costs  of  implanting  flat  and  spherical  geometries  which  correspond  to  the  best 


and  worst  case  geometries,  respectively  It  is  important  to  note  that  these  costs  are  projected  for  a 
state-of-the  art  facility  purposely  designed  for  high  though-put  production  runs  The  costs  of 
64 cents! cm1  for  spheres  and  I4cenrs/cm:  for  flats  are  realistic  estimates  for  this  facility  Much  of  the 
implantation  work  currently  conducted  is  of  an  experimental  nature  or  involves  a  limited  number  of 
work  pieces  The  vendor  is  continually  faced  with  new  processing  problems  involving  the  fabrication  of 
jigs  and  masks  for  new  geometries  This  necessitates  long  set  up  times  which  would  not  occur  for  large 
scale  production  runs.  Therefore,  the  cost  of  implanting  work  pieces  for  research  or  limited  numbers 
of  samples  is  presently  3  to  5  times  those  quoted  in  the  Table.  Accordingly,  vendors  are  presently 
charging  $300  to  $500  per  hour  for  time  that  the  beam  is  on  the  target 


Table  4  6  —  Estimate  of  the  operating  costs  for  the  Man  Tech 
ion  implantation  facility  (from  Smidt  and  Sariwell,  1985) 


1  shift 

2  shifts 

3  shifts 

Capital  investment11 

S100K 

S100K 

S100K 

Full  time  operator 

S  70K 

$175Kd 

$280Kd 

Management 

S30K 

S60K 

$90  K 

l  tilities,  consumables 

S  25K 

$  50K 

$  75 K 

Maintenance  support 

S  20K 

S  30K 

$  40K 

Total  annual  cost 

S245K 

S4ISK 

$58  5  K 

Hourly  cost0 
Throughput  costs'-' 

SI 36/h 

SI  15/h 

S 1 08/ h 

Balls 

$0.64/cm; 

$0. 54/cm, 
$0  1 2/cm2 

SO.  50/cm 

Flats 

$0. 14/cm2 

SO  I  I/cm 

ji  SI  VI/ 10  years  -  SlOOK/year 

hi  2000  h  s^nt l  year  ai  0  9  available  “  1800  h 

ci  2  *  10  lons/cm  .  10  mA 

d1  Assume  lime  and  one  half  for  shift  work 

To  estimate  the  time  for  a  given  sample  area  we  note  that  it  required  16  seconds  to  implant  a  dose 
nt  III1  ions  lor  a  beam  current  of  1  mA.  If  the  total  beam  current  of  the  vendor  is  known,  one  can 
find  the  time  to  implant  the  required  dose  from. 


1 1  <s)  = 


1 6  (mA  -  s ) 
10' 7  (ions) 


(f>  (ions/ cm2 
/  (mA  ) 


x  A  (cm2 


(4 .35) 


where  <!>  is  ihe  required  dose.  I  is  the  total  beam  current  available  at  the  vendor's  facility,  and  A  is  the 
total  area  to  be  implanted.  Following  Smidt  and  Sariwell  this  number  must  be  divided  by  two  beam 
utilization  factors  F,  and  Fh  which  have  values  between  0  and  I.  If  the  total  beam  current  available 
must  be  reduced  to  permit  adequate  workpiece  cooling,  then  the  cooling  factor  /,  will  be  less  than  I. 
Because  in  most  implantations  the  beam  is  raster  scanned  across  a  beam  defining  aperture,  the  beam 
spends  a  fraction  of  time  on  the  apertures  and  masks  rather  than  on  the  sample.  This  factor  is  highly 
dependent  on  the  workpiece  geometry  To  ensure  dose  uniformity  in  the  implanted  area,  the  beam 
must  sweep  entirely  off  the  sample  on  each  side  of  the  aperture  A  typical  beam  spot  diameter  is  1/2 
cm  f  or  an  area  of  I  cm  x  1  cm.  the  beam  is  incident  on  the  mask  for  as  much  time  as  it  is  on  the 
sample  and  />,  =  0  5  For  a  5  cm  *  5  cm  swept  area.  Fh  =  0.83.  Beam  utilization  for  rotating  cylindr¬ 
ical  geometry  lor  a  single  workpiece  is  similar  to  plane  geometries.  However,  if  several  workpieces  are 
to  be  implanted  simultaneously  and  if  masks  are  required  to  limit  sputtering  at  large  incident  angles 
then  a  packing  factor  must  be  introduced  into  Fh  to  account  for  the  mask  areas  For  spherical 
geometries  Fh  is  on  the  order  of  0.1  to  0.2.  Including  the  beam  utilization  Factor  in  the  equation  for 
implantation  time  it  we  have. 
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This  section  on  ion  beam  processing  has  covered  range  statistics,  sputtering,  heat  flow  calcula¬ 
tions,  and  time  and  cost  of  implantation  processing  and  has  included  enough  detail  to  allow  the  reader 
to  plan  the  processing  parameters  for  any  given  ion-implantation  treatment.  As  a  further  aid,  the  next 
section  presents  case  studies  which  use  the  formulae  and  tables  presented  here  in  order  to  illustrate  the 
actual  calculation  of  implantation  parameters. 

V.  DETERMINATION  OF  IMPLANT  PARAMETERS  AND  QUALITY  CONTROL 

Let  us  suppose  that  implantation  of  P  ions  into  flat  steel  specimens  and  implantation  of  Ti  ions 
into  cylindrically  shaped  steel  specimens  is  required  for  a  particular  application.  What  is  the  ion  dose 
required?  What  power  density  limit  is  appropriate  to  avoid  overheating  the  samples'.’  and  so  on.  This 
section  gives  a  step-by-step  procedure  for  the  calculation  of  first-order  estimates  of  ion  dose,  ion 
energy,  beam  power,  time  and  cost  of  implantation,  and  in  so  doing,  illustrates  the  use  of  the  formulae 
given  in  the  previous  section. 

A.  Flat  Geometry  Workpiece 

The  problem  at  hand  is  to  implant  P  ions  into  a  stainless  steel  in  planer  geometry.  The  workpiece 
is  5cm  x  5cm  =  25cm2  in  area  and  is  1  cm  thick.  The  vendor  has  a  maximum  beam  energy  of  90  keV 
and  a  total  beam  current  of  1  mA  available.  It  is  desirable  to  create  an  amorphous  surface  layer  which 
requires  a  concentration  of  20  to  40  atomic  %  of  P.  We  do  not  want  to  overdose  the  sample  because, 
as  we  have  seen,  too  much  dose  can  roughen  the  surface  and  because  too  much  ion  dose  unnecessarily 
raises  the  cost  of  the  implant.  Also,  since  the  material  is  440C  steel,  the  temperature  during  implanta¬ 
tion  cannot  exceed  150°C. 

The  procedure  is  to  calculate  Rp ,  \Rp  S,  <t>s  (90%),  the  beam  power  required  to  maintain  tem¬ 
perature  below  150°C,  and  the  time  and  cost  of  the  implantation  treatment. 


Interpolation  of  Table  (4.1)  for  P  ions  incident  on  an  Fe  target  gives. 


Rp  -  [55-29]  x  Q.8  +  29  =  50  nm 


\R„  -  [23-14]  x  0.8  +  14  =  21  nm. 


From  equations  (4.7-4.10)  we  compute  Sn  =  76  eV.  Substituting  entries  for  Fe  of  U„  =  4.29  eV  and 
N  =0.0848  A  !  from  Table  (4.2),  and  M2  =  56  and  =  31  into  equation  (4.6)  gives, 

S  =  0.042  x  76  x  [0.3 1 9]/ [4.29  x  0.0848]  -  2.8. 


(90%),  C\  (90%) 


Equation  (4.25)  for  the  dose  to  reach  90%  of  sputter  saturation  at  normal  incidence  gives. 


<*>,  (90%)  =  0.96  x  (8.48  x  1022  atm  I  cm2)  x  [50  +  >/2  21](10-?  cm/nm)/ 2.8 

-  2.3  x  1 01 7/  cm7. 


Using  equation  (4.21)  for  the  surface  concentration  at  90%  of  sputter  saturation  yields. 


Cs  (90%)  =  0.9x  1/2.8  =  0.32. 

This  concentration  is  in  the  range  of  20  to  40%  P  that  is  required  so  a  dose  of  2.3  x  10 l7/>/cm2  is 
appropriate.  The  concentration  of  P  is  then  approximated  by  a  32%  concentration  from  the  surface  to  a 
depth  Rp  =  50  nm,  and  0  for  depths  greater  than  50  nm. 

Power  Limit 

The  area  that  the  beam  will  scan  is  25  cm2.  The  total  power  available  is  90  keV  x  1  mA  =  90 
Watts.  Because  the  beam  is  scanned  off  the  samples  on  all  sides  to  achieve  uniformity,  the  beam 
spends  some  time  off  the  samples  so  that  the  true  average  power  is  less  than  90  watts.  It  is  also  neces¬ 
sary  to  determine  the  method  of  heat  sinking  the  samples.  We  will  assume  that  they  are  fastened  with 
a  metal  powder  bearing  glue  that  has  good  heat  conduction  properties  (h=l  W/cm2- K).  Assume  that 
the  entire  available  beam  power  is  incident  on  the  sample  or  90W/25 cm2  —  3.6  W /cm2.  Using  equa¬ 
tion  (4.34)  the  temperature  rise  as  a  result  of  the  interface  between  the  sample  and  the  heat  sink  is 

\  T  =  J/ h  =  3.6  (Wlcm2)l\{W/cm2-K)  =  3.6  K. 

This  is  a  very  acceptable  temperature  drop  across  the  interface  but  notice  that  if  the  samples  were 
only  lightly  clamped  to  the  heat  sink,  the  value  of  h  could  be  <  0.01  and  A  T  could  be  360K.  This 
shows  that  the  method  of  heat  sinking  the  sample  is  extremely  important  for  maintenance  of  tempera¬ 
ture  control.  Even  if  the  samples  are  well  heat  sunk,  they  may  still  be  overheated  because  of  heat  con¬ 
ductance  limitations  through  the  sample  to  the  heat  sink.  In  our  case  the  surface  temperature  at  steady 
state  is  given  by  equation  (4.33),  or 

AT  =  90  ( Wtcm 2)  x  l(cm)/0.21(  W/cm—K)  =  429  K 

where  the  thickness  of  the  sample  is  1  cm  and  the  heat  conductance  k  is  obtained  from  Table  (4.4) 
from  the  entry  for  a  tool  steel.  Thus,  even  with  good  heat  sinking  the  surface  reaches  a  temperature  of 
300  +  429  =  729  C.  To  maintain  a  temperature  limit  of  150°C  requires  that  AT  be  less  than  150°C  - 
25°C  (room  temp.)  =  125  K.  This  can  be  achieved  by  reducing  the  beam  current  by  1  mA  x  125 
K/429  K  =  0.3  mA.  Note  that  if  the  sample  thickness  were  only  0.1  cm,  AT  would  be  43  K  which  is 
an  acceptable  temperature  rise  and  the  full  available  beam  power  could  be  utilized. 

Equation  (4.32)  can  be  used  to  calculate  the  temperature  rise  for  the  worst  case  scenario  where 
there  is  no  heat  sinking  with  the  result, 

T  -  T  =  2.3  x  1017  (ions! cm2)  x  90,000  (eV)  x  1.6  x  10“  19  (Joules/eV)/ 3.7(  Ws/cm'-K)  -  895A 


or.  if  T„  =  25  C,  then  Tma,  =  920  C.  This  is  clearly  too  hot  for  most  applications.  Hence,  implanta¬ 
tion  without  heat  sinking  can  be  considered  only  for  doses  about  one  order-of-magnitude  less,  or  about 
3  x  10lh/ cm'  for  steels. 

Time  and  Cost 

The  time  required  to  implant  this  target  is  computed  using  equation  (4.36)  with  the  result. 


1 6  ( mA  —  s )  ^  2.3  x  IQ1  (ions/ cm2) 
10' 7  (ions)  *  1.0  (mA)  x  1.0  x  0.3 


x  25  (cm2)  =  5 1 


I 


where  Fb  =  1.0  is  estimated  for  the  beam  utilization  factor  and  f,  =  0.3  Fh  is  always  10  when  F,  is 
less  than  1  for  a  flat  geometry  workpiece  or  in  general  when  there  is  no  packing  factor  associated  with 
placing  parts  in  the  beam.  If  the  sample  were  1  mm  thick,  instead  of  1  cm,  or  if  it  had  the  thermal 
conductivity  of  Cu,  then  full  beam  power  could  be  used,  F,  would  be  1.0  and  the  implant  would  take  18 
min.  For  an  hourly  charge  of  $300/hour,  these  samples  would  cost  approximately  $300/hour,  for  a  sin¬ 
gle  batch.  For  a  large  number  of  such  implants  the  cost  could  be  as  low  as  50.14  x  25  cm2  =  53.50  per 
flat  (from  Table  4.6). 

B.  Cylindrical  Geometry  Workpiece 

Suppose  we  have  a  3"  long  tool  steel  cylinder  that  is  clamped  along  a  1"  length  at  one  end  such 
that  there  is  good  heat  conduction  at  the  interface  (h  =  1  H//cm2-K).  The  remaining  2"  length  of  the 
1/2"  diameter  cylinder  is  to  be  implanted  with  Ti  ions  at  an  energy  of  100  KeV  as  it  is  rotated  at  a 
speed  of  10  rpm.  The  treatment  requires  the  sputter  saturated  dose.  The  available  beam  current  is  1 
mA  and  the  vendor  can  either  flood  the  entire  diameter  of  the  rod  or  place  a  mask  in  front  of  the  rod 
to  confine  the  beam  to  a  width  equal  to  the  radius  of  the  rod.  This  limits  the  impingement  angles  to  ± 
30°  with  respect  to  the  surface  normal  whereas  in  the  former  geometry,  all  angles  ±  90°  are  allowed  to 
impinge  upon  the  sample. 

S 

From  Table  (4.1),  Rp  and  \Rp  are  40  nm  and  18  nm,  respectively.  From  equation  (4.6)  the  sputtering 
coefficient  is 

0.042  x  149  [0.271}  =  ? 

^  4.92  x  0.0848 

Sartwell  and  Baldwin  (1985)  have  established  that  during  implantation  of  reactive  ions  into  steels,  sur¬ 
face  carburization  causes  the  effective  sputtering  coefficient  to  be  reduced  by  about  a  factor-of-two  in 
vacuum  pressures  greater  than  10-6  Torr.  Therefore,  a  more  realistic  estimate  of  the  sputtering  coeffi¬ 
cient  is  4.7/2  =  2.4. 

- 

We  will  compute  <t>s  and  Cs  for  implantation  with  and  without  a  mask.  For  a  mask  with  a  width 
equal  to  the  radius  of  the  cylinder,  the  beam  is  confined  to  impingement  angles  ±  30°  from  the  surface 
normal.  The  beam  is  incident  at  angles  smaller  than  30°  for  most  of  the  time  that  the  surface  is  rotat¬ 
ing  through  the  ion  beam  so  that  the  effective  angle  is  much  less  than  30°.  Because  the  cosine  function 
does  not  depart  appreciably  from  a  value  of  1  for  angles  <30°,  it  is  a  good  approximation  to  assume 
normal  incidence  for  this  geometry.  Note  in  Fig.  (4.14)  that  only  a  small  error  is  introduced  by  assum¬ 
ing  ti  =  0  and  ti  =  30°  are  equivalent.  Thus,  for  this  mask  configuration  equation  (4.25)  gives, 

<*,  (90%)  =  yy  (8.48  x  1022)  (40  +  V2  18  (lO  7) 

=  2.2  x  10l7/cm2 

and  the  surface  concentration  from  equation  (4.21)  is, 

C\  (90%)  =  0.9x0.99/2.4  =  0.37. 

Therefore,  the  profile  of  Ti  ions  may  be  approximated  by  a  37%  concentration  from  the  surface  to  a 
depth  of  40  nm 

The  same  calculation  may  be  done  for  the  case  of  no  mask  with  the  help  of  a  simplifying  approxi¬ 
mation  that  the  effective  angle  of  implantation  is  30°.  This  angle  can  be  rationalized  by  noting  that  as  a 


point  on  the  diameter  rotates  through  the  beam,  it  spends  one-half  of  the  time  at  angles  greater  than 
30°  and  one-halt  ot  the  time  at  angles  less  30°  Then  equation  (4.29)  for  </>,  ( 90% )  gives, 

<t>,  (90%)  =  (8.48  x  1  O’2]  (cos  3o]w  [40  cos  30  +  s/2  18]  (lO  7) 

-  1.6  x  101 7 cm2 


and  equation  (4  28)  for  C\  gives, 

C\  =  0.9  x  0.95  x  0.787/2.4  =  0.28. 

Therefore,  a  first-order  approximation  of  the  profile  is  a  Ti  concentration  of  28%  extending  to  a  depth 
of  40  cos  (30°)  =  35  nm.  These  values  of  C\  and  Rp  for  the  two  methods  of  masking  compare  well 
with  the  data  in  Fig.  (4.15)  for  Ti  implanted  AISl  52100  steel.  The  trend  of  smaller  concentration  and 
smaller  depth  in  going  from  normal  incidence  with  a  mask  to  no  mask  is  reproduced  and  the  absolute 
values  are  an  adequate  first  order  estimate.  This  agreement  shows  that  with  reasonable  values  of  Rp , 
±Rr,  and  S,  and  if  other  artifacts  are  not  prevalent  such  as  radiation  enhanced  diffusion  and  preferen¬ 
tial  sputtering,  then  reasonable  agreement  between  measured  and  calculated  profiles  may  be  expected. 

Power  Limit 

In  both  methods  of  masking  the  total  power  on  the  sample  will  be  essentially  the  same  since  all 

ihe  beam  power  available  can  be  concentrated  in  the  rod  whether  or  not  it  is  masked.  There  may  be 

small  differences  in  the  power  because  of  the  beam  utilization  factor  Fh  will  differ  for  the  two  cases. 
Equation  <4  33)  can  be  used  to  calculate  the  maximum  power  than  can  be  delivered  to  the  rod  in  order 
not  to  exceed  150°C.  For  our  steel  rod,  k  =  0.2  and  there  is  good  thermal  heat  sinking  at  one  end.  The 
rod  diameter  D  is  1/2"  and  the  length  /  exposed  to  the  beam  is  2".  First  we  find  if  the  area  of  the 
damping  system  is  large  enough  to  conduct  the  heat  from  the  workpiece,  across  the  interface,  and  into 
the  heat  sink  which  is  kept  at  room  temperature.  The  rod  is  clamped  along  the  1"  section  which  has  a 
surface  area  n  D  /  =  itx  ).27cmx2.54  cm  =  10  cm1.  Equation  (4.34)  indicates  that  for  full  beam 
power  incident  on  the  cylinder  A  T  across  the  interface  will  be 

A  T  =  1°°.  *l\0cm\  _  K 
1  H/cm2 

which  is  a  very  acceptable  value. 

The  maximum  temperature  of  the  rod  lip  is  given  by  equation  (4.33)  multiplied  by  1/2  to  account 

lor  the  fact  that  the  power  is  incident  along  the  entire  rod  (Hirvonen,  1986).  The  temperature  rise  for 

full  available  power  incident  on  the  rod  becomes  A  T  =  Jl/2k  but  the  100  Watts  of  power  must  conduct 
through  the  cross-sectional  area  of  the  rod.  nD2/2.  in  order  to  reach  the  heat  sink,  so  the  effective 
power  density  incident  on  the  rod  for  the  purpose  the  heat  conduction  calculation  is  100  Watts/1.27 
cm2  =  79  Watts/cTw2  The  temperature  rise  is  then 

A  T  =  79  x  5. 08/(2  x  0.2)  =  1000A. 

Therefore,  to  maintain  :ne  tip  temperature  below  150°C.  the  available  beam  current  at  100  keV 
must  be  reduced  by  the  factor  I  mAx  ( 1 50-25)/ 1000  =  0.125  mA  or  an  average  power  of  12.5  Watts 
incident  on  the  sample 

Time  and  Cost 

the  beam  utilization  factor,  is  less  than  l  only  when  the  temperature  rise  is  not  a  problem  and 
available  beam  current  limits  the  implantation  lime,  and  is  less  than  1  only  when  the  temperature 


rise  limits  the  beam  current  on  the  sample.  Therefore,  for  this  case  Fb  —  1  and  F,  =  0.125  because 
the  beam  current  had  to  be  reduced  from  the  1  mA  available  to  0.125  mA.  The  time  of  implantation 
given  by  equation  4  36  for  the  two  methods  of  implantation  becomes 


16  (m.4-s)  2.2  x  101  (ions/cm2)  n  x  1.27  x  5.08  (cm2) 

- - -  x  - - -  x  - 


0.125  x  1.0 


'  101’  (ions)  1.0  (m.4)  0.125x1.0 

=  95mm  (for  ±  30°  mask), 

16  1.6  x  10p  7t  x  1.27  x  5.08  . 

i,  =  — rr  x  - — -  x  — —  — - — — =  69  min  (for  no  mask) 

'  101  10  0.125  x  1.0 


where  the  area  to  be  implanted  is  7tD  /  and  /  is  the  2"  exposed  length  of  the  cylinder.  The  cost  for  the 
former  is  $300/hxl.58h  =  $475  and  for  the  latter  the  cost  if  $300/hxl.l5h  =  $345.  For  batch  process¬ 
ing  the  cost  and  time  may  be  considerably  reduced  because  F,  can  become  1.0  with  Fh  as  high  as  0.8. 
For  cylindrical  geometry  and  batch  processing  the  price  could  be  as  low  as  20  cm2  x  $0.2/cwi2  =  $4  per 
rod  The  rotation  speed  of  10  rpm  indicates  that  for  the  case  of  no  mask,  10  rpm  x  60  min  =  690 
revolutions  of  the  rod  during  implantation.  This  is  an  ample  number  of  revolutions  to  ensure  that 
there  is  uniform  coverage  of  the  rod.  The  total  number  of  revolutions  should  be  greater  than  100  to 
achieve  a  uniformity  of  better  than  1%. 


C.  Quality  Control 


In  general,  ion  implantation  does  not  change  the  visual  appearance  of  a  part  such  as  color,  luster, 
shape  or  roughness.  Unlike  parts  that  have  been  coated  with  TiN  by  chemical  vapor  deposition  and 
which  have  a  gold  color,  the  ion-implantation  vendor  returns  a  part  to  the  customer  which  may  look  no 
different  that  it  did  prior  to  treatment.  Understandably  the  customer,  unfamiliar  with  the  process,  may 
have  doubts  as  to  whether  the  correct  dose  was  used  for  the  implantation  or  whether  the  parts  were 
cooled  properly  and  not  overheated.  Whether  or  not  overheating  occurred  can  be  tested  with 
microhardness  measurements  before  and  after  implantation.  Whether  a  part  received  the  correct  dose 
is  a  quality  control  question  that  can  only  be  addressed  by  experimental  measurements  of  implanted- 
atom  profiles  after  processing. 


There  are  several  implanted  atom  profiling  techniques,  some  of  which  were  used  to  generate  the 
profiles  shown  in  Figures  (4.7  -  4.15).  The  most  common  are  Rutherford  Backscattering  Spectroscopy 
(RBS).  Nuclear  Reaction  Profiling  (NRP),  Auger  Electron  Spectroscopy  combined  with  sputter  erosion 
of  the  surface  (AES),  Secondary  Ion  Mass  Spectroscopy  (SIMS),  and  X-ray  Fluorescence  Spectroscopy 
(XRF).  All  of  these  techniques  are  available  at  various  vendors  that  specialize  in  surface  analysis  ser¬ 
vices.  Table  (5.1)  lists  features  of  each  of  the  methods  and  thus  provides  guidelines  for  which  method 
may  apply  to  a  specific  ion-target  combination. 


Table  5.1  —  Methods  Used  to  Obtain  Implanted  Atom 
Depth  Concentration  Profiles 


IONS 

Detected 

Sensitivity 

Accuracy  of 
Dose  Meas. 

Depth 

Profile  Meas. 

Typical 

Lateral 

Resolution  (mm) 

RBS 

V/,  >  V/: 

good 

good 

yes 

1 

Vf,  <  .vu 

poor 

fair 

yes 

1 

NRP 

V/,<  52 

good 

good 

yes 

1 

AES 

V/,  >  5 

fair 

fair 

yes 

0.1 

SIMS 

V/,  >  1 

excellent 

fair 

yes 

0.1 

XRF 

V,>  24 

good 

fair 

no 

0.01  (in  SEM) 

"  V.' 


RBS  is  overall  the  most  convenient  and  precise  method  when  the  mass  Ml  of  the  ion  is  heavier  than 
the  A/i  of  the  substrate.  NRP  is  the  most  precise  for  light  elements  in  heavy  substrate  but  it  is  more 
time  consuming  than  AES  or  SIMS.  Therefore,  the  precision  of  NRP  is  frequently  used  to  calibrate 
AES  and  SIMS  methods  for  a  given  system  so  that  the  more  convenient  methods  can  be  used  to  gen¬ 
erate  profiles  in  a  shorter  space  of  time.  XRF  is  useful  for  quality  control  during  batch  processing.  If, 
using  one  of  the  other  four  techniques,  it  is  known  that  the  desired  profile  has  been  produced  in  a  part, 
then  XRF  can  be  calibrated  with  this  part  using  x-ray  heights  of  an  element  in  the  substrate  and  the  sig¬ 
nal  from  the  implanted  ion.  This  is  a  measure  of  the  retained  dose  Then  other  parts  treated  the  same 
way  in  a  batch  process  can  be  checked  by  XRF  to  see  if  the  same  dose  is  retained.  Since  XRF  is  a  very 
common  and  easily  used  instrument,  this  is  a  quick  and  inexpensive  quality  control  procedure. 

Another  aspect  of  quality  control  is  the  uniformity  of  ion  dose  over  the  implanted  area.  All  of  the 
methods  in  Table  (5.1)  have  lateral  resolution  capabilities  and  could  be  used  to  profile  different  regions 
of  the  implanted  area  to  check  for  uniformity  of  ion  dose. 

VI.  SUMMARY 

This  chapter  has  touched  upon  all  of  the  practical  issues  involves  in  applying  ion  implantation  to 
parts  and  tools.  In  Section  II,  the  reader  can  determine  whether  or  not  ion  implantation  is  of  interest 
by  applying  the  criteria  and  features  listed  in  Table  (2.3)  and  by  perusing  the  list  of  current  applications 
in  Table  (2.2).  Once  satisfied  that  ion  implantation  is  applicable,  the  fundamental  physical  principles  of 
the  process  are  briefly  described  in  Section  III  for  those  interested  in  the  formation  of  the  surface  alloy. 
This  section  may  be  skipped  with  no  consequence  to  the  understanding  of  the  subsequent  sections. 
Section  IV,  the  heart  of  the  chapter,  is  a  detailed  account  of  the  considerations  required  to  successfully 
ion  implant  a  part  and  includes  tables  and  formulas  for  the  computation  of  the  implanted  atom  concen¬ 
tration  profile,  computation  of  sample  heating  and  methods  to  control  it,  and  computation  of  the 
implantation  time  and  cost.  Section  V  is  a  summary  of  Section  IV  which  illustrates  the  planning  of 
implantation  treatments  in  a  step-by-step  procedure.  Materials  scientists  and  engineers  who  were 
unfamiliar  with  ion  implantation  technology  should  after  reading  this  chapter,  feel  comfortable  discuss¬ 
ing  the  details  of  the  treatment  with  vendors  or  with  researchers  in  the  Field. 

The  primary  advantages  driving  the  development  of  ion  implantation  technology  for  metals  appli¬ 
cations  are  the  low  processing  temperature  and  the  fact  that  the  surface  finish  is  not  effected  (no  refin- 
ishing  required).  It  is  likely  that  surface  coating  technologies  will  never  overcome  the  surface  finish 
advantage  of  ion  implantation,  but  the  low  temperature  processing  advantage  could  conceivably  be 
equaled  by  new  and  innovative  coating  technologies.  For  example,  if  adherent  and  bulk  property  TiN 
CVD  coatings  could  be  applied  at  125°C,  some  of  the  applications  now  using  ion  implantation  may  be 
also  treated  using  this  method.  Another  new  technique  under  development  is  ion  beam  assisted  deposi¬ 
tion  (IBAD)  where  a  low  energy  beam  (100  to  1000  eV)  of  nitrogen  ions  strikes  a  surface  on  which  Ti 
is  simultaneously  being  deposited  by  means  of  a  physical  vapor  deposition  technique.  This  method 
shows  promise  of  being  able  to  produce  bulk  property  refractory  films  at  very  low  temperatures.  In 
spite  of  these  possible  new  developments,  ion  implantation  processing  of  metals  will  continue  to  play  a 
role  in  applications  and  in  materials  research  and  development  because  of  the  speciality  of  the  process. 
By  virtue  of  the  unique  ability  to  implant  any  ion  into  practically  any  target,  ion  implantation  as  a  speci¬ 
alty  industry  can  be  expected  to  grow  in  pace  with  materials  technology  in  general. 
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Fig.  2.1  —  Schematic  view  of  the  development  of  the  concentration  profile  of  ions 

implanted  to  high  doses. 
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Schematic  view  of  ion  implantation  process  (top)  and  a  single  collision 
cascade  (bottom)  (adapted  from  Picraux,  1984). 
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KINETICS  OF  SOLID-SOLID  TRANSFORMATION 


NUCLEATION 

CONTROLLED 

DIFFUSION 

/CONTROLLED 

TIME  FOR  TRANSFORMATION  (LOG  SCALE) 


LOG  TIME  — 

Time-Transformation-Temperature  diagram  demonstrating  effect  of  the  fast-quench 
(10I4K/s)  in  a  collision  cascade  on  formation  of  metastable  microstructures. 
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I  —  The  projected  range,  Rp,  and  standard  deviation  of  the  projected  range 
A/?.,  of  nitrogen  ions  implanted  into  iron  versus  the  initial  ion  energy. 
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Fig.  4.4b  (Continued)  —  Development  of  implanted  ion  distribution  with  increasing  dose  for  cases,  a) 
R*lbRp  *  6  x/2  and  b)  Rp/hRp  «  V2.  W  is  the  amount  of  material  sputtered  away  relative  to  the 
projected  range,  R„  (from  Grabowski  et.  al.,  1984). 
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Fig.  4.5  —  Dose  to  achieve  90%  and  99%  of  saturation  dose  versus  the  modified  ratio  of  Rp  and  A Rp. 
The  curve  labeled  /3  corrects  the  surface  concentration  for  the  tail  of  the  Gaussian  distribution  which 
extends  beyond  the  surface. 
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C  AND  Ti  PROFILES  IN 
TMMPLANTED  STEEL 
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Fig.  4.10  —  Concentration  vs  depth  profiles  for  Ti  and  excess  carbon  in  Ti-implanted 
AISI  52100  steel  at  several  doses  and  energies  (from  Singer,  1984). 
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Fig.  4.13  —  Concentration  vs  depth  profiles  for  Be-implanted  pure  Ni  for  several  ion  doses  and 
energies.  Profiles  were  obtained  with  the  Be  (p,a)  Li  nuclear  reaction  technique  (from 
Westmoreland  and  Hubler,  1977). 
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Fig.  4.21  —  Dimensionless  surface  temperature  rise,  6S,  vs  dimensionless  time  (or  dose),  t,  for 
incident  power  density,  J,  assuming  one  dimensional  heat  transfer.  /  represents  plate  thickness,  or 
cylinder  or  sphere  radius;  k,  thermal  conductivity;  a,  thermal  diffusivity;  4><  ion  dose;  and  E,  ion 
energy.  Heat  is  withdrawn  from  the  target  in  region  III  by  a  heat  sink  (from  Grabowski  and  Kant, 
1983). 
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Fig.  4.23  —  Temperature  rise  of  tool  steel  rod  versus  the  rod  diameter  and  length  for  an  incident  power 
density  of  1  W/cm1  along  the  entire  rod  length.  The  rod  is  assumed  to  be  heat  sunk  at  one  end  over  a 
length  equal  to  one  diameter  with  an  interface  conductance  of  1  (  W/cm2  -  K)  (after  Hirvonen,  1985). 


Fig.  4.24  —  Photograph  of  apparatus  used  to  rotate  and  water-cool  samples 
in  an  ion  implantation  vacuum  chamber  (from  Grabowski  et.  al.,  198S). 


Fig.  4.25  -  Photograph  of  Navy  Manufacturing  Technology  Implantation  Facility  located  at  Spire 
Corporation  From  left  to  right;  Ion  Implanter,  defocussing  section,  work  piece  manipulation 
chamber. 


